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DESIGN  MANUAL 
for 

Self-Acting  Gas  Lubricated  Bearings 
for  Use  in  Electrically  Drivn  Blowers 


I .  INTRODUCTION 

This  manual  is  a  direct  outcome  of  the- requirements 
set  forth  in  BuShips  Contract  NObsr-87522  and  is  intended 
for  use  in  the  design  of  unit  type  self-acting  gas  bearings 
applicable  to  electrically  driven  blowers  in  the  10  to  500 
cfm  range.  The  work  of  preparing  the  manual  was  accomplished 
by  the  Rotron  Manufacturing  Company  with  Mr.  D.  D.  Fuller  serv¬ 
ing  as  Consultant  on  gas  bearings. 

The  information  contained  in  the  manual  has  been  assem¬ 
bled  utilizing  the  existing  available  literature  for  self¬ 
acting  gas  lubricated  bearings,  insofar  as  possible,  the 
data  has  been  verified  by  additional  test  data  accumulated 
at  Rotron  Manufacturing  Company. 

In  its  simplist  form  the  construction  of  a  radial, self¬ 
acting  gas  bearing  can  be  thought  of  as  a  plain  shaft  running 
in  a  very  close  fitting  bore  and  appears  as  a  relatively  simple 
problem. However,  the  clearances  between  shaft  and  bore  must 
be  carefully  established  so  as  to  insure  that  the  two  elements 
are  separated  by  a  gas  film  during  operation.  This  film  carries 
the  applied  load  and  also  acts  as  a  spring  with  a  definite 
spring  rate.  Therefore,  besides  determining  the  load  carrying 
capacity,  the  problem  resolves  itself  into  the  determination 
of  the  spring  rate  and  the  stability  characteristics  of  the 
spring.  As  we  consider  the  parameters  responsible  for  the 
bearing  characteristics,  such  as  the  pressure  distribution  in 
the  gas  film,  the  viscosity  of  the  gas,  and  the  speed  and 
dimensional  considerations,  the  complexity  of  the  problem  be¬ 
comes  more  apparent. 

Before  reviewing  the  general  concepts  of  gas  lubricated 
bearings,  it  would  do  well  to  define  some  of  the  instabilities 
and  terminology  used  in  gas  bearing  design. 

Self-Acting  Bearing:  A  non-pressurized  bearing  that  oper¬ 
ates  on  a  self  induced  gas  film.  It  is  generally  tensed  "hydro- 
dynamic  bearing”  to  distinguish  it  from  the  pressurised  "hydro¬ 
static  bearing". 
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Attitude-Eccentricity  Locust  The  path  that  the  center 
of  the  rotating  component  follows  under  different  load  and 
speed  conditions. 

Synchronous  Whirl t  The  orbital  motion  of  the  center 
of  the  rotating  journal* bearing  combinations  around  the 
center  of  the  stationary  member  in  the  direction  of  rotation 
at  a  frequency  equal  to  that  of  the  spin  speed.  If  the  amount 
of  unbalance  present  is  small,  it  is  possible  to  pass  through 
this  natural  frequency  by  increasing  the  operating  speed  of 
the  unit.  At  this  point  the  spinning  component  will  shift  or 
invert  from  its  geometrical  center  to  its  mass  center  or  more 
specifically  to  its  center  of  gravity.  The  speed  at  which 
this  occurs  is  called  the  "inversion  point"  critical  speed  or 
alternately  the  "synchronous  resonant  whirl”  critical  speed. 

Half-Frequency  Whirl t  The  orbital  motion  of  the  center 
of  the  spinning  member  around  the  stationary  member  in  the 
direction  of  rotation  at  a  frequency  of  one-half  or  less  of 
the  spinning  or  rotational  frequency. 

As  will  become  apparent  in  the  design  section  of  the 
manual,  dimensional  control  of  the  clearances  in  a  bearing 
and  the  type  of  loads  the  bearing  must  sustain  are  generally 
the  key  to  successful  operation  of  the  bearing  once  the  bear¬ 
ing  design  is  established.  This  points  out  the  desirability 
to  design  a  specific  piece  of  rotating  equipment  around  the 
gas  bearing  in  order  to  maintain  the  desired  control  within 
the  bearing.  In  the  case  of  electrically  driven  air  moving 
devices,  bearing  location  within  the  motor,  electrical  unbal¬ 
ance  loads,  air  impeller  location  and  resultant  loads  should 
be  optimized  in  regard  to  their  effect  upon  bearing  perform¬ 
ance  to  insure  reliability  of  the  equipment. 

Use  of  the  Gas  Bearing  Design  Manual 

Insofar  as  possible,  all  curves  and  formulations  con¬ 
tained  in  this  manual  have  been  checked  experimentally  either 
at  Rotron  or  by  the  experiments  conducted  by  some  of  the  refer¬ 
enced  sources  contained  at  the  end  of  this  manual.  It  should 
be  understood,  however,  that  experiments  are  conducted  under 
ideal  conditions  with  known  loads  generally  applied  through  the 
bearing  support  center.  Some  variation  in  calculated  results 
can  be  expected  from  the  need  to  approximate  the  value  of  act¬ 
ual  loads  as  well  as  from  the  variation  in  the  true  clearance 
in  the  bearing  caused  by  taper  or  out-of-round  conditions.  It 
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is  the  attempt  of  this  manual  to  provide  the  designer 
with  the  latest  state-of-the  art  information  on  gas  bear¬ 
ing  design.  Some  good  judgement  on  the  part  of  the  designer 
is  required  when  considering  his  specific  application.  If 
these  considerations  are  understood  the  manual  should  be 
a  very  useful  guide  for  initiation  of  a  gas  lubricated 
blower  design.  Verification  of  the  design  will  necessarily 
be  accomplished  by  actual  test. 

The  recommended  procedure  for  utilizing  the  manual  is 
to  follow  the  order  of  presentation  of  the  manual,  i.e.t 
determine  bearing  size  and  load  capabilities,  and  then  deter¬ 
mine  bearing  stability  characteristics.  An  example  is  given 
for  each  computation  to  assist  the  designer  in  making  use  of 
the  information.  Included  in  the  manual  is  a  section  on 
motor  magnetic  load  influences  of  a  two  pole  alternating 
current  motor  to  assist  the  designer  in  predicting  actual 
loads . 


Review  of  Basic  Concepts  in  Gas  Lubricated  Bearings 

In  order  to  make  effective  use  of  this  design  manual 
it  is  necessary  for  the  designer  to  understand  fully  the 
basic  concepts  of  gas  lubricated  bearings.  Initially,  it 
can  be  stated  that  a  gas  lubricated  bearing  is  a  type  of 
bearing  that  has  its  running  surfaces  completely  separated 
by  a  film  of  gas.  Even  when  load  1b  applied  to  the  bearing, 
the  surfaces  are  kept  out  of  contact  by  this  gas-film  cushion. 

As  with  liquid- lubricated,  fluid-film  bearings,  the  gas 
film  is  established  by  either  of  two  basic  actions.  The 
first  is  characteristic  of  self-contained  bearings  operating 
in  an  open  bath  of  liquid  or  gas.  The  relative  motion  of 
the  sliding  parts  themselves,  generates  a  film  pressure  with¬ 
in  the  clearance  spaces  of  the  bearing  and  insures  separation 
of  the  surfaces  and  load-carrying  capacity.  This  type  of 
bearing  is  called  the  hydrodynamic  (liquid)  or  self-acting 
(gas)  bearing. 

Figure  (1)  is  an  interesting  drawing  taken  from 
Drescher  (1)  and  shows  schematically  this  phenomenon  of 
self-acting  pressure  generation.  An  actual  pressure  trace 
for  an  experimental  journal  bearing  is  shown  in  Figure  (2) 
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taken  from  Ford,  Harris  and  Pantall  (2).  Since  the  vis¬ 
cosity  of  a  gas  is  relatively  low  at  room  temperature, 
about  1/50  of  that  of  water  or  1/3500  of  that  of  an  SAE 
10  petroleum  oil,  the  load-carrying  capacity  resulting 
from  this  self-generated  pressure  is  likewise  relatively 
low.  If  the  bearing  has  to  start  from  rest  under  a  static 
load  and  lift  itself  off  the  "ground0  so  to  speak,  and  be¬ 
come  gas  borne,  unit  loads  are  generally  limited  to  1  to 
2  psi  based  on  the  projected  area  of  the  bearing,  if  the 
load  is  applied  after  rotation  has  been  established  the 
unit  load-carrying  capacity  may  reach  a  practical  maximum 
of  about  lOpsi,  again  based  on  the  projected  area  of  the 
bearing,  although  higher  values  have  been  reported  in  labor¬ 
atory  test  data.  Cole  and  Kerr  (3)  list  average  pressures 
as  high  as  26.5  psi. 

With  the  second  basic  type  of  bearing,  the  externally- 
pressurized  or  hydrostatic  bearing,  the  limitation  on  load¬ 
carrying  capacity  no  longer  exists.  In  this  type  of  bearing 
high  pressure  gas  is  supplied  to  the  clearance  space  in  the 
bearing  and  the  surfaces  are  forced  apart  as  a  result. 

Limits  on  load-carrying  capacity  are  dependent  therefore 
only  on  the  practical  limits  of  gas  pressure  and  gas  volume 
that  are  conveniently  available.  Externally-pressurized 
gas  bearings  have  been  built  with  unit  load-carrying  capa¬ 
cities  as  high  as  50  psi.  There  is,  however,  no  reason 
why  this  figure  of  50  psi  cannot  be  exceeded  if  necessary. 
Since  these  bearings  are  not  dependent  on  speed  for  the 
generation  of  pressure  they  can  conveniently  establish  sepa¬ 
ration  of  surfaces  and  carry  rated  load  with  no  rotation  or 
sliding.  There  is  always, of  course,  flow  of  gas  through  the 
bearing  from  the  high  pressure  source  to  the  lower  anfoient 
pressure.  One  of  the  conditions  for  successful  operation 
of  the  externally-pressurized  bearing  therefore,  requires 
that  this  volume  of  gas  must  be  constantly  supplied,  in 
some  applications  where  relatively  high  static  loads  are 
involved,  hydrostatic  gas  may  be  needed  only  for  starting. 
After  normal  operating  conditions  have  been  established  the 
bearing  may  then  continue  to  operate  as  a  self-acting  bearing. 

However,  in  this  manual  we  will  consider  only  the 
self-acting  gas- lubricated  bearing. 

Most  of  the  complications  that  arise  in  the  analysis 
and  application  of  gas- lubricated  bearings  stem  from  two 
factors,  first,  the  compressibility  of  the  gas  and,  second. 
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Fig.  1  -  Self-Acting  Pressure  Generation  in  a  Journal  Bearing.  (Ref.  1 ). 


the  low  damping  characteristics  of  the  gas  film.  Some  of 
the  bearing  characteristics  that  result  and  a  number  of 
the  problems  associated  with  these  parameters  will  now  be 
discussed. 


When  Harrison  published  his  hydrodynamic  analysis  of 
a  gas- lubricated  journal  bearing  in  1913,(4),  including  the 
effect  of  compressibility,  his  equations  showed  the  effect 
of  ambient  pressure  and  led  to  the  establishment  of  a  dimen¬ 
sionless  group  of  variables  that  seemed  to  be  significant 
in  describing  the  performance  characteristics  of  gas  bearings. 

For  a  journal  bearing  this  wast 

&£§  or  (£)1  (1) 

P.c*  *a  c 
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ANGLE  FROM  LOAD  LINE  (DEGREES) 


Fig.  2  -  Experimental  Pressure  Distribution  in  a  Journal  Bearing.  4x2  Inch  Diameter 
Bearing  Operated  By  a  Radial  Clearance,  C— 0.00 12  Inch  in  Air  At  Atmospheric 
Pressure,  (Ford,  Harris  &  Pantall,  Ref.  2) 


With  the  usual  units,  v  is  the  absolute  viscosity  in 
reynes  (lb.sec/in.2) ;  u  is  the  surface  velocity  of  the  jour¬ 
nal  (in/sec);  r  is  the  radius  of  the  journal  (in.);  Pa  is 
the  absolute  ambient  pressure  (psia) ;  c  is  the  radial  clear¬ 
ance  (in.);  and  u  is  the  rotational  velocity  of  the  journal 
(radains/sec. ) .  See  Figure  (3). 

This  dimensionless  group  of  variables  is  designated  by 
the  Greek  letter  A  ,  (capital  lamda) ,  and  is  termed  the 
"compressibility  bearing  parameter".  The  same  symbol  A  is 
used  for  similar  compressibility  bearing  parameters  for 
other  types  of  gas- lubricated  bearings. 


6 


I 


Fig.  3  -  Journal  in  360°  Full  Bearing. 


For  a  thrust  bearing  the  corresponding  dimensionless 
group  of  variables  isi 


(2) 


Where  1  is  the  length  of  the  pad  in  the  direction  of 
motion  (in.),  h.  is  the  minimum  film  thickness  (in).  Other 
symbols  as  defined  above.  Please  see  Figure  4. 


by* 


For  a  spiral-grooved  thrust  plate,  A  may  be  described 


A 


6uUJ 


h2P, 
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where  j  is  ths  radial  dimension  of  tha  plate  (in);  U  the 
mean  surface  speed  (in/sec);  and  h  the  film  thickness  (in) 
which  in  this  case  is  a  constant.  See  Figure  (5). 


A 

////?/////// 

-  M  __ 

Fig.  4  •  Schematic  Diagram  of  Slider  Thrust  Bearing. 
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Fig.  5  •  Grooved  Thrust  Bearings.  (Ref.  2). 


Zt  has  been  shown  by  Ford,  Harris  and  Pantall,(2), 
by  Ausaan  , (5),  and  others  that  these  dimensionless  para¬ 
meters  designated  by  A  are  pertinent  gas  bearing  parameters 
which  define  the  state  of  operation  of  gas- lubricated  bear¬ 
ings  in  a  similar  fashion  to  the  Soomerfeld  variable  S  for 
bearings  operating  with  liquids  or  other  incompressible 
lubricants. 
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N'  is  the  journal  rotational  speed,  (revs/sec);  Pa  is 
the  average  pressure  based  on  the  applied  load  W,  divided 
by  the  projected  area  of  the  journal,  (W/dl) ,  (psi) ;  d  is 
the  journal  diameter  (in);  other  symbols  as  before. 

The  load-carrying  capacity  of  these  bearings,  as  well 
as  other  characteristics,  can  be  expressed  as  a  function  of 
A  .  For  example.  Figure  (6)  shows  a  curve  for  a  journal 
bearing  at  a  fixed  value  of  eccentricity  ratio  c,  indicat¬ 


ing  the  general  influence  of 
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variable  ambient  pressure. 

When  load  is  applied  to 
a  journal  bearing  the  shaft 
moves  to  an  eccentric  posi¬ 
tion  within  the  bearing 
clearances.  Figure  (3).  The 
eccentricity  ratio  is  a  mea¬ 
sure  of  how  much  the  center 
of  the  shaft  O'  displaces 
itself  from  the  center  of 
the  bearing  0.  Eccentri¬ 
city  ratio  is  defined  as 
the  eccentric  distance  e 
divided  by  the  radial  clear¬ 
ance  c.  or. 


At  light  loads  and  high  speeds,  with  small  unbalance, 
the  eccentric  distance  is  very  close  to  zero  or  e  approaches 
0  and  consequently  c  approaches  0  and  the  shaft  is  nearly 
concentric  with  the  bearing.  At  heavier  loads  and  lower 
speeds  the  shaft  center  becomes  more  displaced  and  when 
metal-to-metal  contact  occurs  between  the  shaft  and  its 
bearing  the  eccentric  distance  e  equals  the  radial  clear¬ 
ance  c  and  the  eccentricity  ratio  c  is  then  1.  The  limits 
on  the  magnitude  of  eccentricity  ratio  are  therefore  be¬ 
tween  zero  and  1. 
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From  a  consideration  of  the  variables  included  in  A , 
one  can  gain  an  understanding  of  the  pertinent  parameters 
that  influence  the  performance  of  self-acting,  gas-lubri¬ 
cated  bearings.  The  ambient  pressure  enters  the  analysis 
since  the  bearing  acts  as  a  pump  and  like  all  pumps,  the 
ratio  of  its  discharge  pressure  to  supply  pressure  is 
limited.  Thus  a  reduced  ambient  pressure  and  density 
around  the  bearing  will  result  in  a  lower  maximum  pressure 
in  the  film.  Naturally  this  means  less  load-carrying  cap¬ 
acity. 


Speed  is  also  a  factor  in  gas- lubricated  bearings  as 
it  is  in  liquid- lubricated  bearings  but  there  is  one  dis- 
tict  difference.  In  the  liquid- lubricated  bearing  the  load 
capacity  continues  to  rise  directly  as  the  speed  increases, 
with  no  limit,  assuming  constant  viscosity  of  the  lubricant. 
With  a  gas  bearing,  however,  although  the  load  capacity  does 
increase  with  speed,  the  rate  of  increase  falls  off  at  high 
speed  due  to  compressibility  effects,  and  the  load-carrying 
capacity  approaches  a  limiting  value.  At  extreme  speeds 
the  load-carrying  capacity  becomes  independent  of  viscosity 
and  speed  and  is  only  a  function  of  the  ambient  pressure 
surrounding  the  bearing.  Figure  (7)  is  from  Pantall  and 
Robinson,  taken  from  (6). 

In  general,  dimen¬ 
sionless  plots  of  load¬ 
carrying  parameter  W' 
against  A  for  various 
values  of  eccentricity 
ratio  may  be  shown  as 
in  Figure  (8),  taken 
from  (2).  Design  curves 
of  this  type  are  in¬ 
cluded  in  the  design  sec¬ 
tion  of  this  maiual. 

SPflO 

Fig.  7  -  Load  Carrying  Capacity  of  Self-Acting  Gas-Lubricated  Journal  Bearing 
As  Influenced  By  Speed.  (From  Pantall  and  Robinson,  Ref.  6) 
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Fig.  8  -  Generalized  Experimental  Correlation  of  Dimensionless  Load  Parameter 
W'vs  A  .  for  1  /d=  2  .  Departure  from  Straight  Line  Relationship  When  A  >1 
Showing  Compressibility  Effect.  (Ford,  Harris  &  Pantall,  Ref  2). 


It  will  be  noticed  that  for  A  greater  than  one,  the 
plot  no  longer  yields  a  straight  line  relationship.  The 
deviation  from  the  straight  line  indicates  the  point  at 
which  compressibility  effects  begin  to  exert  their  influence 
It  is  generally  agreed,  that  for  journal  bearings,  the  com¬ 
pressibility  and  variable  density  influences  are  negligible 
as  far  as  load-carrying  capacity  is  concerned  for  values  of 
A  somewhat  below  1  and  the  classical  theory  for  incompress¬ 
ible  lubricants  may  be  used,  if  desired,  without  incurring 
a  significant  error. 

This  does  not  apply  to  attitude  angle  however,  (angle  * 
of  Pig. 3),  and  this  will  be  discussed  in  the  design  section. 

Another  point  must  be  made  regarding  ambient  pressures 
in  the  very  low  range,  approaching  those  of  a  vacuum.  The 
present  theory  of  gas- lubricated  bearings  assumes  that  the 
gas  can  be  considered  as  a  continuous  fluid.  This  assump¬ 
tion  is  no  longer  valid  at  very  low  pressures  where  the 
molecular  mean  free  path  of  the  gas  becomes  comparable  in 
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magnitude  to  the  film  thickness  in  the  bearing.  When 
this  happens ,  effects  occur  which  can  no  longer  be  ex¬ 
plained  by  continuum  flow  theory.  Slip  occurs  at  the 
boundary  between  the  bearing  surface  and  the  gas.  There¬ 
fore,  instead  of  assuming  vanishing  velocities  at  the 
walls,  slip  velocities  must  be  introduced.  This  has  been 
done  by  Burgdorfer, (7) . 

If  the  ratio  of  minimum  film  thickness  to  mean  free 
path  is  called  n,  it  can  be  shown  that  for  values  of  n 
less  than  100  a  noticeable  effect  on  bearing  pressure  and 
load- carrying  capacity  may  be  expected.  A  typical  graph 
is  shown  in  Figure  (9)  from  (7).  This  graph  happens  to  be 
for  a  slider  bearing.  A  similar  condition  would  exist  for 
a  journal  bearing. 

Typical  values  of  mean  free  path  at  atmospheric  pres¬ 
sure,  aret 

Hydrogen  4.43  x  10” ^  inches 

Helium  7.32  x  10~^  inches 

Air  2.52  x  10“*?  inches 

Neon  5.20  x  10~°  inches 

Thus  for  n  *  100  for  air  at  atmospheric  pressure,  the 

film  thickness  at  the  entrance  to  a  tilting  pad  bearing 
would  be  of  the  order  of  2.52  x  10~4  inches  (0.000252  inch). 

Such  values  are  not  uncommon  for  self-acting  bearings, 
even  at  atmospheric  pressure.  For  low  ambient  pressures 
the  effect  would  be  even  more  marked  and  should  be  considered 
in  any  analysis  and  design  for  these  operating  conditions. 

This  can  be  illustrated  by  a  table  (Table  I)  taken  from  (7) 
where  the  ratio  of  radial  clearance  to  molecular  mean  free 
path  was  computed  for  a  number  of  cases  as  reported  by  wild- 
mann,  (8).  It  can  be  seen  because  of  either  the  very  small 
radial  clearances  in  these  journal  bearings  or  because  of 
the  low  ambient  pressure,  the  value  of  this  ratio  in  all  cases 
is  less  than  100. 

Another  and  different  series  of  problems  is  associated 
with  gas- lubricated  bearings  because  of  the  low  damping  pro¬ 
perties  of  the  gas  films.  Because  of  this  characteristic, 
certain  dynamic  instabilities  need  to  be  anticipated  and  con¬ 
trolled  with  machinery  operating  on  gas- lubricated  bearings. 
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Fig.  9  -  Pressure  Distribution  for  Compressible  Fluids  Along  a  Straight 
Slider  Bearing.  (Ref.  7) 

With  liquid- lubricated  bearings  these  might  have  been 
suppressed  or  they  might  have  passed  by  unnoticed  be¬ 
cause  of  the  greater  damping  action  of  liquids.  With 
gas- lubricated  bearings,  however,  these  possible  insta¬ 
bilities  must  be  investigated,  and  in  most  cases  design 
modifications  can  be  made  to  avoid  them  or  at  least  to 
mitigate  their  action. 

The  various  types  of  these  dynamic  behavior  character¬ 
istics  will  be  described  and  analysed  in  the  design  section 
of  this  report.  The  most  important  instabilities  are  half¬ 
frequency  whirl  and  synchronous  whirl. 
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II.  DESIGN  OF  GAS- LUBRICATED  JOURNAL  BEARINGS 


A.  Load-carrying  Capacity  with  Static.  Uni-Directional  LoadB 

A  typical  journal  bearing  is  shown  in  Figure  (3).  This 
is  a  schematic  view  with  a  greatly  exaggerated  clearance, 
running  at  some  angular  speed  under  the  action  of  a  load  W. 
When  starting,  with  solid  contact  between  the  shaft  and  the 
bearing  at  a  point  directly  under  the  load,  the  shaft  will 
start  to  run  up  the  right  side  of  the  bearing.  After  rising 
a  few  degrees  of  arc  it  will  enter  a  region  where  the  gas 
comes  between  the  bearing  and  the  journal.  Slip  will  begin. 
Then  as  the  journal  comes  up  to  speed  it  will  revolve  faster 
and  faster  and  build  up  a  wedge-shaped  film,  which  will  then 
place  the  journal  in  the  position  shown  in  Figure  (3). 

A  polar  diagram  of  film  pressures  is  shown  in  Figure  (3) 
representing  the  positive  film  pressures  in  the  loaded  zone 
of  the  bearing.  There  will  be  negative  pressures  (below 
ambient  pressure)  in  the  unloaded  zone  of  the  bearing  which 
will  also  contribute  to  the  load-carrying  capacity.  These  are 
not  shown  in  Figure  (3). 

In  general,  the  pressure  generation  in  the  clearance  of 
a  journal  bearing  is  described  for  liquid  lubricants  by  the 
Reynolds'  equation.  Ref. (9), 

h  +  16  ohl  (6> 


where  p  is  the  pressure  at  any  point  in  the  film 

h  is  the  film  thickness  at  any  circumferential 
location 

v  is  the  absolute  viscosity 
U  is  the  surface  speed  of  the  journal 
x  is  the  coordinate  in  the  circumferential 
direction 

Z  is  the  coordinate  in  the  axial  direction, please 
see  Figure  (10) 

For  a  gas,  the  compressibility  has  to  be  included.  This  is 
represented  by  the  mass  density  being  a  variable  rather  than 
a  constant  as  with  the  liquid  case  so  that  the  characteris¬ 
tic  equation  becomes.  Refs.  (10)  and  (11), 
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(7) 


General  Solution  for  Load  Capacity 

Exact  solutions  for  Equations  (6)  and  (7)  cannot  be 
obtained  so  that  approximate  methods  have  been  used  to  devel¬ 
op  design  information.  The  best  of  these  solutions  have  re¬ 
sulted  from  the  application  of  digital  computer  techniques. 

The  work  of  Elrod  and  Malanoski  (12)  and  (13),  Elrod  and  Burg- 
dorfer  (10),  Raimondi  (14),  Hays  (15)  and  Gross  (11)  will  now 
be  presented  as  a  general  design  approach  to  360°  journal  bear¬ 
ings. 

Figures  (11)  through  (15)  show  the  load-carrying  para¬ 
meter,  w/dlpa,  plotted  against  the  compressibility  bearing 
parameter,  A,  for  various  values  of  eccentricity  ratio,  c. 

W  is  the  applied  load,  pounds 
d  is  the  journal  diameter,  inches 
1  is  the  journal  length,  inches 
Pa  is  the  ambient  pressure,  psia 
A  is  defined  by  Equation  (1) 

«  is  defined  by  Equation  (5) 
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On  each  sheet  is  included  the  table  of  data  from  which  the 
particular  figure  was  drawn.  There  is  one  figure  for  each 
1/d  ratio. 
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LOAD  PARAMETER 
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Fig.  1 1  -  Theoretical  Load-Carrying  Parameter  vs  Compressibility  Parameter 
for  Full  Journal  Bearing  for l/d=co(Ref.  10  and  14). 
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TABLE  II 


Compressibility  Parameter 
vs 

Load  Parameter  for  1/d  of  Infinity 
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TABLE  III 


Cotnprsssibility  Paranstsr 
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KARIM*  COMPRESSIBILITY  PARAMETER  A* 

FOR 
(RKP.  14) 

Fig.  13  •  Theocertcd  Load-Carrying  Parameter  vj  Compressibility  Parameter  for 
Full  Journal  Bearing  for  l/d=2  (Ref.  14) 
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BEARING  COMPRESSIBILITY  PARAMETER  A« 
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(RIF  14) 

Fig.  1 4  •  Theoretical  Load-Carrying  Parameter  vs  Compreuibility  Parameter  for 
Full  Journal  Bearing  for  l/d=l  (Ref.  14) 
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pm  pm  om  om 


TABLE  V 


Compreaaiblity  Parameter 

va 


A 

Load  Parameter 

for  1/d 

-  1 

W/dlPa 

W/dlPa 

c 

A 

0 

0 

0.4 

0 

0 

0.24 

0.00902 

0.24 

0.0413 

0.6 

0.0223 

0.6 

0.101 

1.2 

0.0429 

1.2 

0.192 

3.0 

0.0862 

3.0 

0.386 

6.0 

0.116 

6.0 

0.545 

12.0 

0.131 

12.0 

0.643 

- 

0.1594 

m 

0.7966 

0 

0 

0.6 

0 

0 

0.24 

0.0186 

0.24 

0.0770 

0.6 

0.0457 

0.6 

0.188 

1.2 

0.0875 

1.2 

0.355 

3.0 

0.176 

3.0 

0.720 

6.0 

0.238 

6.0 

1.047 

12.0 

0.273 

12.0 

1.280 

" 

0.3335 

m 

1.6349 

0.8 

0 

0 

0.06 

0.0404 

0.24 

0.171 

0.6 

0.428 

1.2 

0.783 

3.0 

1.555 

6.0 

2.276 

12.0 

2.835 

m 

3.6652 
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BEARING  COMPRESSIBILITY  PARAMETER  A-  ——  [rl* 

» *4t-k  Pb  w 

(REF.  14) 

Fig.  1 5  •  Thaoratical  Load-Carrying  Paramatar  vs  CompranibSity  Paramata r 
for  Full  Journal  Baaring  for  l/d=£  (Rif.  14) 
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TABLE  VI 


Compressibility  Parameter 
vs 

Load  Parameter  for  1/d  ■  jj 


A 

*«lpa 

c 

A 

W/dlPa 

0 

0 

0.4 

0 

0 

0.24 

0.00288 

0.24 

0.0141 

0.6 

0.00719 

0.6 

0.0350 

1.2 

0.0143 

1.2 

0.0693 

3.0 

0.0348 

3.0 

0.163 

6.0 

0.0634 

6.0 

0.288 

12.0 

0.0974 

12.0 

0.449 

24.0 

0.122 

24.0 

0.587 

m 

0.1594 

m 

0.7966 

0 

0 

0.6 

0 

0 

0.24 

0.00599 

0.24 

0.0290 

0.6 

0.0149 

0.6 

0.0721 

1.2 

0.0297 

1.2 

0.141 

3.0 

0.0717 

3.0 

0.325 

6.0 

0.129 

6.0 

0.560 

12.0 

0.199 

12.0 

0.873 

24.0 

0.252 

24.0 

1.163 

ee 

0.3335 

ee 

1.6349 

0.8 

0 

0 

0.06 

0.0184 

0.24 

0.0768 

0.6 

0.201 

1.2 

0.391 

3.0 

0.834 

6.0 

1.342 

12.0 

1.991 

24.0 

2.602 

ee 

3.6652 
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Before  proceeding  with  an  example  we  will  need 
some  viecotiy  data  for  air.  It  should  be  understood 
that  knowing  the  viscosity,  the  solutions  presented  in 
this  report  are  applicable  to  other  gases  as  well  as  air. 
Table  VII  is  a  summation  of  data  from  References  (16)  and 
(17). 


TABLE  VII 


Viscosity  Data  for  Air 

Absolute  Viscosity  in  Reyns  (lb. sec/in2 )x  10  9 


Pressure  psia-» 
Tern]?.  °F 

14.7 

294 

735 

1000 

70 

2.65 

2.75 

2.88 

2.95 

122 

2.86 

2.895 

2.962 

3.02 

212 

3.18 

3.205 

3.24 

3.275 

302 

3.465 

3.48 

3.50 

3.52 

Thus  the  viscosity  of  air  at  70°F  and  14.7  psia  is 
2.65  x  10“  reyns. 

Figure  (16)  is  a  plot  of  the  viscosity  data  from  Table  VII. 


Fig.  1 6  -  Absolute  Viscosity  of  Air.  (Refs.  ( 1 6)  and  ( 1 7). 
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LOAD  APPLIED  TO  BEARING  W(LBS) 


uai  capacity  far  Boar laps  with  A  »  X.  mppoao 
w  viiM  to  totoniM  tto  fiUlf  acaaotriclty  rati*  ail  ai»i> 
■mm  fila  tRlckooaa  foe  a  totrtof  with  tto  fallavlap  apaci- 

f icatiaaat 


Rotor  Malpht 
•pood 

IlkiMt  riMMTt 
LaBricatlap  rUii 
Aofciaot  TMparator* 

Be it iap  LaoptR 

Boar lap  Diaaotor 

Boar lap  Radial  cloaraaca 


0.118  18.  (Cor  aoa  Roar  Up) 

24.000  ROM 

14.7  polo 

Air 

100°T 

0.50  la* 

0.50  la. 

0.00025  ia. 


Roariav  Caloalatiaoai 

•-(“if*1 

•  •  l.U  a  10**  nyH  Pl«.  (U) 

251)  raOa/aoc. 

r  -  0.25  iocAoa 
»,  -  14.7  paia 
c  •  0.00025  iacAoa 

*  -  toupy«Aj«u 

•*'  I  •  itMt  • 1000  *  »•* 

, .  om&icijt  mu  « t»* 

1  -  1.71 


0  0J  CXt  0J  0l4  OlS  OlO  0.7 


ECCENTRICITY  RATIO  € 


Fig.  17  *  Theoretical  Load-Carrying  Capacity  360°  Journal  Bearing. 
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For  1/d  •  1,  and  *  •  3.75.  In  Fiquro  (14)  wo  My  bo  oblo  to  ootn> 
blioh  tho  operating  occontrlclty  ratio,  one#  we  have  ovolootod 
*/dl»a. 

V  -  O.Ut  lb. 

dlfa  -  0.5  x  0.5  x  14.7  •  3.675  lb. 

Thon  -  0.110/3.675  *  0.0321 

Thoa  in  Fiyuro  (14)  with  W/dlFa  -  0.0331  and  *  •  2.75  we 
find  that  tho  operating  *  haa  a  valuo  of  somewhere  bo two an 
‘  ■  0.1  and  «  ■  0.  Probably  a  valuo  of  •  •  0.05  lo  about  aa 
cloao  aa  It  can  bo  oatlxtatod. 

Tho  minimum  flla  thlcknoaa  la  thon,  frea  Figure  (3) 

hc  •  e  U-O  (0) 

-  0.00035  I 1-0.05) 

-  0.00025  x  0.05 

-  0.0002375 

h0  •  0.00024  lnehoo 


Thia  bearing  would  need  to  be  checked  for  at  leaat  half-frequency 
whirl  atability  and  for  the  poaaibility  of  aynchronoua  vibration 
before  the  design  would  be  conaidered  aa  complete.  (See  appropri¬ 
ate  sections  in  manual) .  If  the  bearing  were  subjected  to  a  range 
of  steady- state  loads,  as  for  example  if  the  unit  were  being  sub¬ 
jected  to  a  wide  acceleration  spectrum,  then  the  performance  pattern 
could  be  calculated  by  assuming  different  values  of  w/dlPa,  for  the 
same  speed  and  ambient  pressure  conditions.  Thus,  in  this  case 
A  *  a  constant  of  value  2.75 


Assumed  Acceleration  w(lbs) 
Condition.  "G's" _ Resultant 


w/dlPa 


c 


1 

2 

3 

4 

5 

6 
7 

e 

9 

10 

20 

30 

These  results  are  plotted 


0.118 

0.0321 

0.05 

0.236 

0.0642 

0.075 

0.354 

0.0816 

0.1 

0.472 

0.1285 

0.15 

0.590 

0.1605 

0.2 

0.708 

0.1927 

0.24 

0.8255 

0.2245 

0.29 

0.944 

0.257 

0.31 

1.061 

0.289 

0.35 

1.18 

0.321 

0.37 

2.36 

0.642 

0.58 

3.54 

0.963 

0.69 

Figure 

(17). 
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Journal  Bwring  Solution*  for  Low  Values  of  A 


It  can  be  observed  in  Figures  (11)  through  (15)  that 
for  low  values  of  A  ,  such  below  a  value  of  one,  the  prob¬ 
able  accuracy  that  can  be  obtained  from  the  use  of  these 
figures  is  considerably  reduced.  It  has  also  been  shown 
by  Ford,  Harris  and  Pantall  (2),  (Figure  8)  that  for  jour¬ 
nal  bearings  the  compressibility  and  variable  density  in¬ 
fluences  are  negligible  as  far  as  load-carrying  capacity  is 
concerned  for  values  of  A  somewhat  below  one  and  the  class!' 
cal  theory  for  inccaipressible  lubricants  aay  be  used. 

This  has  also  l>een  shown  very  clearly  by  Raiaumdi  (14) 
and  Figure  (18)  is  reproduced  froa  his  paper. 


COUfftCMIttllTY  PARAMKTCII 
*•*£-&)* 


Fig.  18*  Discrepancy  h  Load  Capacity  Between  Theories  Based  On  Incompressible 
and  Compressible  Lubricants.  (Ref.  1 4). 


For  convenience  therefore,  and  in  soae  instances  for 
increased  accuracy,  a  solution  based  on  incompressible  lub¬ 
ricants  aay  be  used.  One  of  the  best  of  these  is  due  to 
Hays,  Ref.  (15).  His  solution  is  based  on  the  use  of  the 
Soaaerfeld  variable,  S,  Equation  (4).  For  the  case  ofA*  0, 
we  have  the  following  relationships  between  c,l/S,  and  1/d 
ratios 
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0.1 

- 

0.01 

0.042 

0.102 

0.228 

0.710 

6.5 

0.2 

- 

0.065 

0.165 

0.398 

0.882 

2.625 

9.5 

0.3 

- 

0.17 

0.363 

0.870 

1.885 

5.305 

14.2 

0.4 

0.15 

0.3 

0.625 

1.485 

3.14 

8.352 

19.0 

0.6 

0.25 

0.55 

1.30 

3.019 

6.08 

14.60 

25.5 

0.8 

0.42 

0.9 

2.09 

4.755 

9.18 

20.31 

31.0 

1.0 

0.65 

1.3 

2.92 

6.513 

12.16 

25.25 

39.0 

1.25 

0.9 

1.75 

3.91 

8.56 

15.46 

30.35 

48.0 

2.50 

1.65 

3.4 

7.25 

15.17 

25.33 

44.0 

59.0 

5.0 

2.25 

4.65 

9.52 

19.49 

31.42 

51.85 

66.0 

10.0 

2.95 

5.5 

10.69 

21.69 

34.50 

55.80 

71.0 

Fran  Equation  (4) , 

1  -  Ii23  (£) 

S  "  h*  * 


Convert ing  H’  -  j?  Pavg  - 


va  have 


1  W/dl  2».c.2 
a  •  M —  (7> 


6  pm 

Multiplying  by  7  and  — *  wa  hava 
6  Pa 


1  _  W/dl  2*  ,c 


a  w/  ua  /v»  w  ¥  « 

S  M  -  V  X  6  p, 

—  ■  (— )  x  — —  x  12* 

S  6W*V  dlPa 

■i  -  4  *  W*  12*  whara  W' 
8  A 
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ECCENTRICITY  RATIO,  € 

Fig.  1 9  -  I/S  Plotted  Against  Eccentricity  Ratio  for  A  Approaching  Zero. 

(Data  from  Ref.  1 5) 


Fig.  20  -  I/S  Plotted  Against  Eccentricity  Ratio  for  A.  Approaching  Zero. 

(Data  from  Ref.  1 5) 


or 


1  12  W'  w 

S  "  A 


(10) 


This  provides  the  relationship  between  the  Sommerfeld  num¬ 
ber  s  and  the  compressibility  parameter  A.  (for  leer  values 
of  A) 


Data  from  Table  VIZI  are  plotted  in  Figure  (19)  and 
the  lower  left  hand  corner  of  Figure  (19)  is  expanded  in 
Figure  (20)  to  permit  more  accurate  evaluation  at  very  low 
values  of  eccentricity  ratio  e. 


jmsatlS-li  ua  ccapute  thi  load-carrying  capacity  of  a 
gaa- lubricated  journal  bearing  for  which  *  la  laaa  than  one.  Air 
at  70  p  and  14.7  pain  la  to  bt  ueed. 

*  ■  2.6s  «  lo"®  froa  Pig.  16 
1  -  0.75  in. 
d  -  0.375  in. 
c  -  0.00025  in. 

■  «  3150  RPM 

«  -  4  ouncaai Pavg  -  1S»0TJ73»U.73  "  P»i 
O  •  61.9  ln/aac. 
r  ■  0.375/2  -  0.1675  in. 

PB  »  14.7  paia 

ccapute  value  of  4  .  Proa  aquation  (1), 


6x2.65xl0-S  61.9  *  0.1675 
(0.00025)'  x  14.7 


*  -  0.201 

Row  the  1/d  ratio  1*  0.75/0.375  -  2. 

On  flrat  figure  (13)  which  include#  tha  affect  of 
preeeiblllty  and  for  a  given  4  and  W/dlfa  find  a. 


dlj£  “  6. 37^3^75x14. 7  "  °'°604 

Than  froa  figure  (13),  a  »  0.4 

the  a)a(aaa  fila  t’  iokneea,  aquation  (6),  la 

h#  »  c  U-»l 

•  0.0002S  (1-0.4) 

•  0.00035  a  0.6 
h0  -  0.00015  ia. 
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Vast  uee  the  naya  eolution  for  inccapreaaible  lubricanta 
correeponding  to  **0.  from  aquation  (») 


1  -  U.37 


Bow  in  rlgura  (It),  with  1/a  *  11.37  and  1/d  »  2,  wa  find  that 
<  •  0.4,  which  ia  the  ana  aa  calculated  above  frcai  riqura  (13). 


Zacidantally,  the  alternate  fora  of  1/a  frea  aquation  (10) 
ylelda 


Therefore 


J  where  W*  •  >010604  from  above, 

j  “  11.33  the  aaae  raault  aa  froa  Bq.  (9). 


Kx ancle  3i  Auaean,  Ref. (18),  citea  experimental  data  frea 
a  nuaber  of  aourcea.  The  firat  la  given  below  froa  acheinberg. 
Data  are  for  an  1/d  ratio  of  1.1  and  *  -  3.5. 


Calculated  calculated 

Calculated 

KxparUwntaX 

c 

c 

a 

w/dl». 

c 

u/a-x) 

(1/4*2) 

By  Interpolation 

0.33 

0.33 

0.4 

0.3 

0.39 

0.385 

0.3 

0.58 

0.44 

0.37 

0.83 

0.63 

0.66 

0.56 

0.63 

1.13 

0.73 

0.73 

0.63 

0.74 

The  agreement  la  quite 

aatiafactory. 

hleo  froa  (18)  we 

have  data  froa 

atemlicht  and  Blwell  for 

•  -  1.3 

and  1/d  -  1.5t 

Calculated  calculated 

Calculated 

Departmental 

t 

c 

f 

B/dl», 

a 

(1/4-1) 

(l/d-9) 

By  Interpolation 

0.11 

0.13 

0.33 

0.13 

0.3 

0.36 

0.19 

0.47 

0.33 

0.393 

0.4 

o.ss 

0.63 

0.43 

o.sas 

0.3 

0.63 

0.68 

0.S1 

0.595 

0.68 

0.83 

0.76 

0.63 

0.699 
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Ry.  21  -  Moment  Capacity  of  Full  Journal  Bearings.  {Amman,  Ref.  5). 


Sow  rspsst  that*  last  calculation*  of  *  using  rigur*  (19)  and  tha 
approach  baaad  on  inooatpraaaibla  lubricants  avan  though*  is  1.3 
and  is  higher  than  tha  raccaatandad  Halt  of  ona. 


W/dl», 

4 

CBlCulStSd 

ri«.  is 

t 

0.11 

3. IS 

0.14 

0.15 

0.26 

7.S4 

0.32 

0.39 

0.4 

11.6 

0.46 

0.55 

O.S 

14.5 

0.53 

0.65 

0.6S 

IS. 7 

0.63 

0.82 

Agra  want  ia  only  fair,  aapaclally  at  highar  valuaa  of  <  . 

Moment-Carrying .Capacity  of  Journal  Bearing  (Over¬ 
hung  Static  Load) 


To  estimate  the  moment- carrying  capacity  of  a  full 
journal  bearing  we  may  use  the  analysis  of  Ausman  (18). 
There  is  little  experimental  verification  of  these  pre¬ 
dictions.  However,  the  results  may  be  used  to  show  what 
bearing  variables  are  significant  in  producing  a  moment¬ 
carrying  characteristic. 


Figure  (21)  is  reproduced  from  Ausman's  paper.  In  it 
is  plotted  a  dimensionless  parameter  H,  against  another 
dimensionless  parameter  which  we  can  call  K.  This  is  done 
Cor  various  values  of  1/d  ratio  of  the  bearing. 


Fig.  22  -  Schematic  Representation  of 
Journal  Bearing  Sustaining 
An  Applied  Moment,  M. 


H 


K 


ka 


6)mr* 

kapac' 


(11) 


l3  6  (12) 


coefficient  of 
polytropic  expan¬ 
sion  of  gaseous 
lubricant  which  for 
our  purposes  is  iso¬ 
thermal  air  so  that 
ka  -  1.0 


M  ■  moment  applied  to  the  bearing  (in. lbs) 


j  ■  misalignment  angle  or  angular  displace¬ 
ment  of  journal  (see  Fig. 22)  radians 
Other  variables  as  before. 
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E»«nplt  4 1  conaldar  ■  Journal  bearing  with  tha  following 
apecif leatlona i 


1  -  1-1/e  in. 

d  -  0.375  In. 
it  ■  1.63  oa. 

R  -  3150  RPM 
t.  -  14.7  paia 
e  •  0.000375  In. 

a  •  1. 65  x  10~9  (Air  at  14.7  pala  a  70°) 


Subatltuta  In  equation  (11)  for  H, 


H 

R 


6x3,65  a  10-9[~80)x  b 
1.0x14.7  x  0.00037 3 J 


2 


R  -  0.0093 


Prom  Plgura  (21)  wo  uaa  a  llnaar  lntarpolatlon  for  low  valuaa  of  H. 


Row  K  la  glvan  aa  Equation  (12),  aolvlng 
for  N  wa  hava 


°-0267  • 

M  >  73.44 


Thua  tha  nowant  dapanda  upon  tha  pamlaalbla  angular  aiaalign- 
nant  4  batwoan  tha  eantral  axla  of  tha  Journal  and  tha  cantral  axla 
of  tha  baarlng.  Ona  can  look  at  tha  altuatlon  In  althar  of  two  wayai 
by  eonaldarlng  tha  4  that  would  raault  f row  tha  application  of  a 
known  nonant  N,  or  oatlnata  tha  naxlnun  angular  dlaplacanant  that 
night  ba  tolaratad  and  lntarprat  thla  In  tarna  of  a  naxlnun  allcwabla 
nonant  that  night  ba  appllad  to  tha  Journal  baarlng. 


Por  thla  agunpla  lat  ua  uaa  tha  lattar  approach.  For  tha 
eoncantrlc  poaltlon  of  tha  Journal  In  Ita  baarlng  tha  flln  thick- 
naaa  at  point  I,  Plgura  (22), will  ba  aqual  to  tha  radial  claaranca 
c,  or  0.000175  inchaa.  Wa  night  than  aaauna  a  doflactlon  aa  ahown 
with  a  nlnlnun  flln  thlcknaaa  at  point  1,  Plgura  (22),  of  0.00005  In. 
and  a  naxlnun  flln  thlcknaaa  at  point  C  of  0.00070  Inchaa. 

Tha  total  langth  of  tha  baarlng  la  1-1/0  In.  ao  that  tha  angla 
to  bo  uaad  baconaa 


tan  4 


0.0007-0.00005 

1.125 


tan  4  •  0.000570  or  vary  cloaaly, 
4  •  0.000570  radlana 


Row  aubatltutlng  la  tha  aquation, 

N  -  73.44 

N  -  73.4  X  0.000570 
N  •  0.0424  Inch  lb  a. 
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B.  Evaluation  of  Friction  in  Full  Journal  Bearings 


Of  interest  in  some  designs  is  the  power  loss  in  the 
bearing  due  to  friction  resulting  from  the  viscous  drag 
effects  of  the  gas.  In  gas  bearings  the  friction  is  de¬ 
scribed  from  the  friction  factor  cj  which  is  defined  as 
follows  t  Mj 


11 


_•  _  2tuUr  1 
cf  “  c - 


(13) 


where  Mj  is  the  friction  moment  in  inch  pounds,  other  sym¬ 
bols  the  same  as  before. 


The  friction  moment  being  considered  arises  from  the 
hydrodynamic  film  forces  acting  on  the  journal  of  a  bear¬ 
ing  during  steady  state  operation,  designated  M,.  As  evi¬ 
denced  on  Figure  (23) ,  when  a  load  WA  is  applied  through 
the  center  of  the  bearing  0,  the  reaction  load  w  (which  is 
equal  and  opposite  to  WA)  occurs  at  the  journal  center  O' 
creating  a  moment  equal  to  WOO'  sin  +  .  If  the  lubricant 
film  in  the  bearing  clearance  space  is  considered  a  free 
body,  then  the  friction  moment  of  the  journal  acting  on 
the  lubricant  film  will  be  and  the  friction  moment  of 
the  bearing  will  be  Mj,.  Summing  up  the  clockwise  and 
counterclockwise  moments  we  havet 


Mjj  +  WOO'  sin  *  -  Mj 


from  Equation  (5)  00' 


Fig.  23  *  Relationship  Between 
Bearing  Friction  Torqve  and 
Journal  Friction  Torqve  (9). 


>  «  c  T  Mj  -  Mjj+Wtc  sin  ♦ 

It  should  be  noted  that  the 
journal  friction  is  greater  than 
the  bearing  friction  if  friction 
measurements  are  being  made  through 
the  bearing.  Of  consideration  also 
is  the  fact  that  the  journal  fric¬ 
tion  varies  with  the  angle  4,wfcich 
is  called  the  attitude  angle  of  the 
bearing  or  simply  attitude,  as  trail 
as  with  the  eccentricity  ratio  c. 
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Fig.  24  -  Friction  Factor  for  Full  Bearing  of  Infinite  Length  (Ref.  1 5). 


The  friction  loss  of  s  360°  journal  bearing  is  rather 
independent  of  A  at  values  of  «  up  to  about  0.5.  This 
swans  that  for  low  c  the  values  of  friction  as  obtained  for 
the  case  of  incompressible  lubricants  (A  ■*  0)  may  be  used. 
This  can  be  seen  graphically  for  exasple  in  Figure  (24) 
which  is  reproduced  from  Ref.  (12)  for  the  infinitely  long 
bearing . 

Table  IX  is  a  compilation  of  data  for  the  infinitely 
long  bearing  for  all  values  of  A  up  to  8.  Values  for  A  *  0 
were  obtained  from  References  (15),  (14),  and  (9). 

TABLE  IX 


Friction  Factors  for  Full  Bearings  of  Infinite  Length 


0.0  0  1.0 

0.2  0  1.09 

0.4  0  1.33 

0.6  0  1.82 

0.8  0  2.87 

0.9  0  4.30 

0.2  0.9622  1.0518 

0.4  0.8209  1.2346 

0.6  0.5418  1.6536 

0.8  0.2052  2.6722 

0.9  0.0715  4.0505 

0.2  1.979  1.0331 

0.4  1.859  1.1528 

0.6  1.460  1.4501 

0.8  0.6348  2.3268 

0.9  0.2221  3.6303 

0.2  4.012  1.0242 

0.4  3.976  1.1090 

0.6  3.643  1.3090 

0.8  2.379  1.8987 

0.9  1.017  2.9412 

0.2  8.058  1.0214 

0.4  8.110  1.0956 

0.6  7.781  1.2648 

0.8  6.256  1.7220 

0.9  4.246  2.4733 
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For  bearings  of  finite  length  a  correction  must  be 
applied  to  the  value  of  cj  ae  obtained  from  Table  IX  and 
Figure  (24).  The  value  of  the  correction,  although  cal¬ 
culated  for  the  case  where  A  ■  0,  may  also  be  applied  to 
bearings  with  values  of  A  up  to  about  6.  A  graph  of  such 
correction  factors  is  shown  in  Figure  (25).  The  data  for 
this  graph  were  obtained  from  Ref.  (15).  Cf  is  the  friction 
factor  for  the  journal  bearing  of  finite  length. 

Raimondi  (14)  also  shows  the  influence  of  compressi¬ 
bility  effects  on  the  friction  of  finite  journal  bearings. 

In  summary,  the  following  rules  may  be  used  to  obtain 
values  of  the  friction  factor* 

1.  For  0  ■  A  ■  6  use  Elrod's  values  of  c^.  Fig.  24, 
corrected  for  1/d  ratio  from  Hays'  Figure  25. 

2.  For  6  ■  A  «  12  use  c^  ■  cj  from  Fig.  24.  No  cor¬ 
rection  for  1/d  ratio  being  required. 

3.  For  A  ■  12  use  c_  «  (Ref.  13). 

f  1-c2 


Kx—ple  5 i  Suppose  we  consider  the  beering  of  Exastple  1. 
A  -  2.75,  1/d  -  1,*  -  0.05 


From  Figure  (24),  0}  ■  1.0 

cf 

Frost  Figure  (25),  the  correction  factor  3  *  1.0.  No  correction  et 

‘  low  values  of  * . 


Therefore,  c{  •  1.0  x  cj 

Cf  -  1.0 

Then  frost  Equation  (13) 

c*  "  2»sor4l/c 
2»MPr2lcs 

or  c 


(14) 


For  this  exaaple,  v  ■  2.68  x  10  reyns 

0  ■  2510  rads/sec  x  0.25  in. lb."  627.5  in/sec. 
r  *  0.25  in. 

1  -  0.50  in. 
c  -  0.00025  in. 


Substituting  in  Equation  (14), 


„  .  **  3.6telQ-9x6a7.5at(0.25)2x0.50»a.O 
"j  0.00025 

2t«  0.260x0. 6275x0. 0625x0.50 
"j  25 

Mj ■  0.00132  inch  lb. .friction  torquo. 

Or  translatod  into  horsepower  at  a  apoad  of  24.000  KM, 


f.h.k. 


33,000 


2skO.00132ii24.000 
12  x  33,000 


It  is  always  possiblo  to  gst  a  good,  ordsr-of ■  magnitude 
chock  on  ths  friction  in  a  journal  boaring  by  using  Eewton's 
definition  of  absolute  viscosity.  As  shown  in  Roforsneo  (9) 
Equation  (6),  page  8, 

M  -  JITp  (15) 


K  ■  absolute  viscosity  in  reyns 
(lbs.  sec/in2) 

A  ■  the  swept  area,  2*rl(in2) 

U  ■  the  surface  speed  of  the  journal 
(in/sec) 

h  ■  the  radial  film  thickness  which  is  c 
Thus  using  the  parameters  from  Example  5, 


N  ■  0.00132  inch  lb.  friction  torque 

In  this  case  the  answer  is  identical  to  the  previous  result  ob¬ 
tained  for  Example  5  since  the  only  difference  between  Equations 
15  and  16  is  the  absence  of  the  c<  factor  from  Equation  (15).  Thus 
when  Cf  •»!,  both  equations  yield  the  same  result. 
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C.  Attitude  Angle  in  Full  Gas- Lubricated  Journal  Bearings 


It  is  clear  from  the  discussions  outlined  m  the  intro¬ 
duction  that  the  journal  center  does  not  coincide  with  the 
bearing  center  under  steady  operating  conditions.  There  is 
actually  a  locus  of  points  that  the  journal  center  will  fol¬ 
low  depending  upon  the  bearing  geometry.  Figure  (26)  shows 
a  typical  locus  of  motion  of  the  center  of  a  journal  in  the 
clearance  of  a  self-acting  journal  bearing.  The  angle  i  is 
the  attitude  angle  formed  between  the  direction  of  applied 
load  and  the  line  joining  the  bearing  and  shaft  centers. 

This  angle  is  shown  in  Figures  (23),  (26)  and  (27)  and  is  of 
importance  in  establishing  both  the  mathematical  and  physi¬ 


cal  criteria  for  stability  in 
journal  bearings. 


w 


Fig.  26  -  Typical  Attitude  -  Eccentricity 
Locus  for  the  Motion  of  the  Center  of  a 
Journal  in  the  Clearance  of 
o  Self-Acting  Bearing. 


hydrodynamic  (self- acting ) 


Attitude  Relationships  for  Journal  Bearing. 


To  determine  the  parameters  affecting  attitude  angle 
it  is  necessary  to  examine  the  hydrodynamic  pressure  profile 
in  the  bearing.  In  Figure  (27)  an  applied  load  is  sub¬ 
jected  to  the  bearing  through  the  bearing  center  0,  and  the 
reactive  load  is  taken  at  the  journal  center  O'.  The  result¬ 
ant  film  force  W  has  a  component  along  the  line  of  centers 
00*  equal  to  W  cos  4. 
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The  hydrodynamic  pressure  at  any  point  in  the  film 
is  a  variable,  p,  a  function  of  speed,  viscosity,  clearance, 
angular  position  in  the  clearance,  eccentricity  ratio,  etc. 

W  cos  4  is  then  opposed  by  that  component  of  the  pres¬ 
sure  film  force  acting  along  the  line  of  centers  00'. 

0-6' 

w  cos  4  ■  f  p  cos  ©  r  1  d©  (16) 

0-0 


And  in  the  same  way  for  the  component  of  the  result- 
and  film  force  W  perpendicular  to  the  line  of  centers  00', 
we  have, 

©-©• 

W  sin  4  -  f  p  sin  ©  r  1  d©  (17) 

0-0 

For  purposes  of  this  explanation  sideflow  or  varia¬ 
tion  of  p  in  the  axial  direction  is  neglected. 


Dividing  Equation  (17)  by  Equation  (16)  yields. 


Wsins 

Wcos4 


©-©' 


i 


p  sin  ©rid© 


©-©' 


(f-0 


p  cos  ©rid© 


tan  4 


(18) 


or  in  general,  tan  4  is  a  function  only  of  eccentricity 
ratio  c  .  This  function  is  unique  for  any  particular  bear¬ 
ing.  For  example,  with  a  liquid  lubricated  journal  bearing, 
assuming  a  pressure  arc  of  180°,  (Ref.  19), 


tan  4 


w  Vl-c* 

2  c 


(19) 


Evaluating  Equation  (19)  we  have  Table  X. 
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TABLE  X 


Attitudg  Angles  for  Liquid  Lubricated.  180°  Journal  Bra. 


c 

tan  4 

4  degrees 

0 

oo 

90 

0.1 

15.65 

86.35 

0.2 

7.70 

82.6 

0.3 

5.0 

78.7 

0.4 

3.6 

74.5 

0.5 

2.72 

69.8 

0.6 

2.095 

64.5 

0.8 

1.18 

49.7 

0.9 

0.762 

37.3 

0.95 

0.517 

27.35 

0.98 

0.319 

17.7 

Whan  plotted  on  polar  coordinate*  a*  shown  in  Figure 
(28) ,  the  resultant  curve  is  known  as  an  attitude-eccen¬ 
tricity  locus  and  represents  the  path  along  which  the  center 
of  the  journal  will  always  traverse. 

In  a  similar  fashion,  the 
attitude  angles  for  certain  gas 
lubricated  bearings  have  been 
calculated  including  the  compressi¬ 
bility  effects  of  the  gas.  Thus 
the  attitude  angles  for  gas- lubri¬ 
cated  bearings  become  a  function 
of  A  as  well  as  c.  Figure  (29) 
shows  attitude  angle  patterns  for 
full,  gas- lubricated  journal  bear¬ 
ings  as  found  in  Ref.  (14).  Notice 
the  effect  of  A  on  these  locii. 
These  diagrams  are  similar  to  Fig¬ 
ure  (28)  where  the  eccentricity 
ratio, c  -  1,  is  replaced  by  the 
radial  clearance. 

For  quantitative  evaluation  of  these  angles,  charts  from  Elrod, 
Ref.  (13),  are  used  since  they  cover  a  greater  range  of  A 
values  with  higher  accuracy  and  also  include  one  additional 
1/d  ratio.  Tables  of  these  values  are  also  included. 


Fig.  28  -  Attitude-Eccentricity  Locus 
for  Liquid-Lubricated, 

1 80*  Journal  Bearing. 
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Fig.  29  -  Attitude-Eccentricity  Locus  Diagrams 
for  Full,  Gas-Lubricated  Journal  Bearings.  (Ref.  1 4). 


lag  having  an  1/d  retie  of 
1  end  e  velue  of  ■  ■  0.6, 
the  following  ettitude 
englee  ere  obtained  from 
Figure  (33)  or  Table  XIV i 


0 

*0' 

0.1 

85 

o.a 

79 

0.4 

69 

0.6 

62 

1.0 

52 

4.0 

25 

8.0 

15 

20.0 

60.0 

3.9 

100.00 

2.7 

The  dependency  ef  etti¬ 
tude  angle  ♦  on  *  ia 
obvious. 


It  can  be  eessn  in  these  figures  that  considerable  dif¬ 
ference  exists  in  the  attitude  angle  between  fluid  theory 
for  incoonpressible  lubricants  (A*0)  and  fluid  theory  for  con- 
pressible  lubricants,  even  for  very  small  values  of  the  coat 
pressibility  parameter  (or  bearing  number) A .  This  is  not  like 
the  case  for  load-carrying  capacity  where  there  is  little  dif¬ 
ference  in  predicted  values  between  liquids  and  gaseous  lubri¬ 
cants  up  to  a  value  of  A  about  one. 

For  gas-lubricated  bearings  the  value  of  the  attitude 
angle  4  must  always  be  determined  through  Figures  (30)  to  (34) 
by  means  of  the  compressibility  parameter. A. 
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1 


0 


0.1  OS  0.4  04  I  t  4  t  •  to 
BEARING  COMPRESSIBILITY  PARAMETER  A* 

OATA  FROM  TABU  XI 


tO  40  00  00100 

T  If 


Rg.  30  •  Attitude  Angle  vs  Compressibility  Parameter  for  l/d=oo  (Ref.  13). 
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TABLE  XII 


Attitude  Angle,#  (deg.)  for  1/d  »  » 


A 

c-*0 

e-0.2 

0 

90 

90 

0.1 

85.58 

- 

0.24 

79.50 

- 

0.5 

68.95 

67.58 

1 

52.76 

- 

2 

34.35 

33.43 

6 

14.75 

14.20 

10 

9.87 

9.42 

20 

5.82 

5.86 

24 

5.09 

5.08 

30 

4.32 

4.29 

40 

3.52 

3.46 

50 

3.01 

2.94 

60 

2.65 

2.58 

70 

2.39 

2.32 

80 

2.18 

2.11 

90 

2.02 

1.95 

100 

1.88 

1.82 

• 

0 

0 

0.4 

0.6 

0.8 

90 

90 

90 

64.13 

55.80 

38.43 

31.87 

27.43 

19.40 

13.45 

11.72 

8.77 

8.80 

7.57 

5.65 

5.32 

4.49 

3.28 

4.59 

3.84 

2.78 

3.84 

3.19 

2.28 

3.08 

2.52 

1.77 

2.60 

2.10 

1.46 

2.27 

1.82 

1.25 

2.03 

1.62 

1.10 

1.84 

1.46 

0.98 

1.69 

1.34 

0.89 

1.57 

1.23 

0.82 

0 

0 

0 

52 


oi  as  04 


4  i  in 


40  oo  ooioo 


■CARING  COMPRESSIBILITY  PARAMETER  A* 
MT*  PROM  TMU  XX 


Pig.  31  •  Attitude  Angle  v*  Compressibility  Parameter  for  1/d-IT  (Ref.  13). 


33 


TABLE  XIII 


Attitude  Angle,  ♦(deg.)  for  1/d  ■  2 


A 

c+0 

c  “0.2 

0 

90 

90 

0.06 

- 

- 

0.1 

86.55 

- 

0.24 

81.77 

81.28 

0.6 

70.17 

68.98 

1.2 

54.46 

53.22 

3 

30.58 

29.94 

6 

18.26 

17.74 

12 

11.00 

11.33 

24 

6.745 

6.715 

30 

5.794 

5.720 

40 

4.780 

4.675 

50 

4.130 

4.014 

60 

3.671 

3.552 

70 

3.327 

3.208 

80 

3.059 

2.940 

90 

2.842 

2.726 

100 

2.662 

2.548 

m 

0 

0 

0.4 

0.6 

0.8 

90 

90 

90 

- 

- 

76.56 

78.80 

71.94 

52.48 

64.47 

54.58 

36.37 

48.63 

39.99 

26.60 

27.47 

22.92 

15.94 

16.10 

13.43 

9.370 

10.20 

8.486 

6.068 

5.907 

4.750 

3.254 

4.998 

3.976 

2.687 

4.050 

3.180 

2.110 

3.456 

2.686 

1.757 

3.042 

2.346 

1.518 

2.737 

2.097 

1.343 

2.501 

1.905 

1.210 

2.311 

1.752 

1.105 

2.155 

1.627 

1.019 

0 

0 

0 

54 


ATTITUDE  ANGLE  ♦  (DEGREES) 
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TABLE  XIV 


Attitud«  Anql«. i  (d»q.)  for  1/d  ■  1 


A 

e+0 

c-^.2 

0 

90 

90 

0.06 

- 

- 

0.1 

88.38 

- 

0.24 

86.11 

85.62 

0.6 

80.35 

79.33 

1.2 

71.24 

69.32 

3 

49.96 

47.93 

6 

31.60 

30.33 

12 

18.81 

17.95 

24 

11.77 

11.66 

30 

10.21 

10.02 

40 

8.528 

8.279 

50 

7.430 

7.168 

60 

6.648 

6.383 

70 

6.058 

5.798 

80 

5.594 

5.340 

90 

5.218 

4.970 

100 

4.904 

4.663 

«• 

0 

0 

0.4 

0.6 

0.8 

90 

90 

90 

- 

- 

77.77 

83.53 

76.92 

54.01 

74.43 

68.09 

37.36 

62.19 

47.96 

27.87 

41.11 

30.50 

17.94 

26.01 

19. 55 

11.90 

15.32 

11.58 

7.140 

9.924 

7.500 

4.687 

8.474 

6.343 

3.909 

6.950 

5.146 

3.118 

5.987 

4.397 

2.630 

5.312 

3.878 

2.297 

4.810 

3.494 

2.053 

4.419 

3.197 

1.865 

4.104 

2.959 

1.716 

3.844 

2.763 

1.594 

0 

0 

0 
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Rg.  33  •  AftHude  Angle  vi  CompranibWy  for  twee  ter  for  l/d*l  (Ref.  1 3). 
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TABLE  XV 


Attitude  Angle.*  (deg.)  for  1/d  ■  h 
A _ c-»0  c  "0.2 _ 0.4 _ 0.6 _ 0.8 


0 

0.06 

0.1 

0.24 

0.6 

1.2 

3 

6 

12 

24 

30 

40 

50 

60 

70 

80 

90 

100 


90 

89.48 

88.75 

86.88 

83.75 
74.79 

61.54 
43.04 
26.03 
21.93 
17.71 

15.15 

13.44 
12.21 
11.26 
10.50 

9.866 

0 


90 


88.54 

86.37 

82.78 

72.45 

58.62 
40.81 

24.62 
20.32 

16.75 
14.47 
12.87 
11.68 
10.74 

9.987 

9.368 

0 


90 


87.48 

83.73 

77.82 

63.63 

48.78 

33.32 

20.24 

16.95 

13.85 

11.91 

10.55 

9.546 

8.760 

8.136 

7.615 

0 


90 


83.41 

74.44 

63.16 
45.53 
33.29 
22.72 
14.07 

12.16 
9.844 
8.405 
7.414 
6.686 
6.119 
5.668 
5.298 

0 


00 

81.51 

61.56 

43.46 

32.60 
22.36 

16.60 

11.74 

7.450 

6.898 

5.507 

4.659 

4.082 

3.661 

3.336 

3.078 

2.868 

0 
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100 


40  40  00.00 


BEARING  COMPRESSIBILITY  PARAMETER  A" 

DATA  FROM  TABU  XV 


Fig.  34  -  Attitude  Angle  vs  Compressibility  Parameter  for  l/d=£  (Ref.  13). 
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In  addition  to  unidirectional  constant  radial 
forces,  the  second  most  common  type  radial  load  that  a 
journal  must  carry  is  mass  unbalance  of  the  rotating  com¬ 
ponent.  Such  an  unbalance  produces  a  radial  force  that 
rotates  at  shaft  speed  as  depicted  on  Figure  (35)  with  a 
magnitude  dependent  upon  both  the  amount  of  unbalanced  mass 
in  the  system  and  the  rotational  speed. 


As  with  a  steady,  unidirectional  load,  a  rotating  force 
will  produce  a  displacement  of  the  journal  within  the  bear¬ 
ing.  In  this  case,  however,  the  center  of  the  journal  will 
not  be  displaced  to  a  fixed  point  in  space  but  will  describe 
a  somewhat  circular  orbit  because  of  the  angular  motion  of 
the  force.  It  is  clear  that  this  angular  motion  occurs  at 
the  same  frequency  as  the  shaft  rotational  speed  and  there¬ 
fore  it  is  generally  termed  ' synchronous  whirl ' .  If  all 
other  radial  forces  were  removed  from  the  journal  then  the 
center  of  the  circular  orbit  of  the  shaft  center  would  corre¬ 
spond  closely  to  the  center  of  the  bearing  as  seen  on  Figure 
(36) path  A.  On  the  other  hand,  if  a  steady  unidirectional 
radial  load  acts  on  the  journal  to  displace  the  journal  cen¬ 
ter  to  some  point  O'  (Fig. 36),  then  the  center  of  the  circu¬ 
lar  orbit  due  to  unbalance  would  also  shift  to  O'  and  the 
new  path  described  by  the  journal  center  would  approximate 
path  B. 


Fig.  35  -  Fore#  F,  Resulting  from  Unbalance  Path  B  Due  to  a  Steady  Load 

of  Rotor  Shaft  Rotating  with  Shaft  Plw*  Unbalooc#- 

At  Synchronous  Speed. 
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In  order  to  determine  the  minimum  bearing  film  thick¬ 
ness  of  a  journal  bearing  containing  unbalance  of  the  rotat¬ 
ing  member  it  is  necessary  to  determine  the  radius  e'  of  the 
circular  orbit  of  the  shaft  center. (See  Fig.  36).  For  small 
values  of  e*  the  eccentricity  ratio  c  and  e'  My  be  calcu¬ 
lated  as  though  the  rotating  load  were  a  stationary  load  of 
the  same  magnitude.  This  has  been  shown  to  be  quite  accurate 
if  c  does  not  exceed  a  value  of  about  0.4  to  0.5.  Thus  the 
analysis  to  predict  the  load-carrying  capacity  of  a  static 
load  may  be  applied  to  n  load  due  to  rotating  unbalance  with 
satisfactory  accuracy  provided  that  the  value  of  eccentricity 
ratio  c  does  not  exceed  about  0.4  or  0.5.  (Ref.  20).  For 
eccentricity  ratios  greater  than  these  values  dynamic  re¬ 
sponses  occur  in  the  gas  film  involving  squeeze  film  forces 
and  other  time  dependent  actions  so  that  the  static  equiva¬ 
lency  described  above  becomes  more  inaccurate. 

To  review  the  accuracy  of  considering  the  unbalance 
load  as  a  fixed  load  the  following  example  taken  from  Ref. 

20  by  Sternlicht  and  Elwell  is  presented. 


ESSBEifi-Z1  *  3*0°  gaa- lubricated,  MU-actiq  Journal  boar¬ 
ing  wea  toatod  on  a  vortical  shaft.  Attitude  angle  and  eccentri¬ 
city  ratio  ware  naaaurod  for  the  abaft  ufeilo  la  operation,  for 
the  following  conditiona,  the  eccentricity  ratio  tor  a  rotating 
unbalance  nay  be  confuted  assming  the  analyaie  for  a  ateady  load 
condition  ia  applicable i 


1/d  -  1 
Dimeter  -  2* 
Length  ■  2* 

Spaed  ■  6000  rpn 


Radial  Clear. ■  0.000904 
Viacoaity  *  2.61xl0",reyne 
Onbal.  S  6000  rpn  ■  12.5  Ib/brg 
Anb. pressure  -  14.7  paia 


A  m  2.63 

_ ,  —8  e  — Hal U  0.2125 

and  14.7x2x2 


Uaing  figure  (14)  or  Table  V,  the  calculated  value  of  i  «  0.27. 

The  neaaured  value  of  c  e  0.275  indicating  excellent  agreanent.  To 
naka  a  noro  conplote  analyeia  eeveral  other  pointa  were  calculated 
with  the  following  resultai 


Speed 

(run) 

W/*.dl  * 

CAlmlatod 

c 

Measured  t 
(Ref  .20) 

2,000 

1.5 

0.0255  0.076 

0.00 

O.OS 

4,000 

5.05 

0.0050  1.755 

0.15 

0.15 

6,000 

12.5 

0.2125  2.63 

0.27 

0.275 

0,000 

22 

0.374  3.505 

0.37 

0.405 

10,000 

35 

0.505  4.30 

0.47 

0.53 

The  reaulte  are  shorn  in  figure  (37)  ccnpared  with  exper 1- 
anntal  data. 


61 


Fig.  37  •  Orbital  Eccentricity  Ratio  for  Journal  Subjected  To  An  Unbalanced 
Radial  Load.  Experimental  Data  from  Ref.  20. 

establishment  of  Limit*  for  Rotating  Unbalance 


The  method  of  analysis  just  discussed  also  pro¬ 
vides  a  means  for  establishing  the  limits  on  unbalance  that 
can  be  tolerated  in  any  particular  rotating  assembly.  This 
limit  will  be  measured  in  terms  of  the  eccentricity  ratio 
that  the  unbalance  produces  and  the  corresponding  minimum 
film  thickness. 


Iwb1«  6i  Aa  an  illuatration  of  thia,  racon alter  baapl*  7 
that  haa  juat  boon  deacrlbed  aa  having  an  accantridty  ratio  of 
«  ■  0.27.  Tha  apaclf lc  boaring  In  taagla  7  had  a  langth  and  dia¬ 
meter  of  2  inchoa  aach  with  a  radial  claaranco  of  0.000504  inchoa. 

Zf  it  la  dealrablo  to  na  la  tain  tha  niniana  filn  thicknaaa  of  thia 
boarlngat  0.00017  lachaa  at  6,000  IN  then  tha  maximon  accantrielty 
ratio  allowable  nay  be  calculated  ualng  aquation  (S) 

h_ 

h0  -  c  [1-tlor  1 - 5  "  * 

-  1  -  .734 

t  -  0.266  or  approx inataly  0.27 

Tor  thia  accantrielty  ratio  and  tha  calculated  bearing  ccm- 
praaaiblllty  number  *  -  2.63,  tha  maxinam  load  nay  be  detenu  load 
frcn  Figure  (14)  to  be  12.3  pounda  par  bearing.  Thia  force  la 
tha  reeult  of  aa  unbalance  which  can  be  determined  free  the  relation 

r unbalance  "  I  * 

where  W  la  the  magnitude  of  the  unbalance, a  la  the  rotational  apaad 
in  radlana  per  aec.,  *g*  ia  the  gravitational  conatant  and  r  ia  the 
effective  radlua  of  the  unbalance. 
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Thus  (or  thia  aaaaqpla,  substituting  tho  uatricil  vsluss,  _  m 

havai 


ia.s  --Al  *  a»l*  *  * 

«r  >  0.0122  lb. in. 
or  Hr  “  0.105  ot.ii.  of  talaUaet 

Tbs  allowablo  aabUiact  Itctu  a  (pact lea  of  boarlao  claar- 
aaco,  alaa  and  apttd,  for  a  given  lubricant  viscosity,  and  la  neat 
caaea  would  bo  establish**  by  the  Halts  oa  alalaaa  fila  ttldattt 
that  weald  bo  acceptable,  la  a  aajority  of  appllcatloaa  a  alalaaa 
fila  tbickasaa,  reealtiag  frca  unbalance .  of  paabapa  free  1/4  to  2/3 
of  the  radial  clearance  would  probably  be  allowable,  i.e.  c  ■  .25 

to  t  •  0.34. 


Synchronous  Rt»en«at  Whirl ( Inversion- Point  critical 

Speed) 

If  the  gas  fila  in  the  bearing  is  considered  ss 
a  spring  supporting  the  rotating  eleaent,  then  any  disturb¬ 
ing  force  will  create  a  deflection  of  this  spring.  In  the 
case  of  an  unbalance  force,  the  disturbing  force  varies 
angularly  indirection  at  a  frequency  equal  to  that  of  the 
rotating  speed  of  the  rotor.  At  any  specific  point  in  the 
gas  fila  the  unbalanced force  thus  acts  as  a  vibrating  force 
with  a  frequency  equal  to  the  rotational  speed.  The  rotat¬ 
ing  shaft  systea  which  is  considered  supported  on  a  gas 
spring,  has  a  natural  frequency  at  which  it  will  resonate. 

If  the  frequency  of  the  unbalanced  vibrational  force  corre¬ 
sponds  to  the  lowest  natural  frequency  of  the  shaft  and  sup¬ 
port  systea,  a  forced  vibration  phenaaenon  occurs  which  is 
frequently  tensed  "synchronous  resonant  whirl”  or  "inversion 
point  critical  speed".  Both  of  the  terns  iaply  that  the 
rotor  is  acting  as  a  rigid  body  and  is  vibrating  on  the  rela¬ 
tively  soft  gas  "springs"  in  the  bearings.  (Ref.  21) 

Extensive  observations  have  been  aade  on  this  type  of 
instability  (Ref.  22).  By  aeans  of  capacitance  probes  Moni¬ 
toring  the  notion  of  the  journal  in  its  bearing,  it  can  be 
shown  that  below  this  resonant  "critical  speed"  the  shaft- 
rotor  systsa  is  rotating  about  its  geoaetric  axis.  After 
passing  through  this  resonant  condition,  the  shaft- rotor  sys- 
tsa  begins  to  spin  about  its  center  of  gravity  axis.  This 
inversion,  or  shift  in  the  center  of  rotation,  can  be  observed 
on  the  screen  of  an  oscilloscope  and  hence  the  origination  of 
the  tera  "inversion-point  critical  speed”. 
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Translator^;  Synchronous  Whirl  Rotating  Journal 


Tha  roaponao  of  a  gaa  film  spring  is  non¬ 
linear,  consequently,  there  is  no  spring  "constant"  but 
rather  a  spring  rate  which  will  increase  with  eccentricity 
ratio  for  this  type  of  bearing.  The  actual  spring  rate 
therefore,  at  a  given  eccentricity  ratio,  may  be  obtained 
by  drawing  a  tangent  to  the  curve  of  load-carrying  capa¬ 
city  plotted  against  eccentricity  ratio  for  a  given  bear¬ 
ing  at  a  fixed  speed  (Refs.  22  and  23) .  Such  a  curve  would 
be  typical  of  Figure  (17). 


rifur*  (1?)  ii  redraw  In  Figure  (38)  including 
m  an  additional  abscissa  tha  alni—  tila  thickneee  eorraaponding 
to  tha  various  valaea  of  * .  Thane  era  calculated  (ran  aquation  (•) 

h#  -  c  [1-4 

For  tha  bearing  daecribed  in  Figure  (17),  e  •  0.00023  inchea. 


QuSOMS  MNNNt  OMM  MMI  QBMM  0 
MINIMUM  FILM  THICKNESS  -  !*•  (Indies) 


Fig.  38  •  Replot  of  Data  from  Example  1  and  Rg.  1 7. 


If  the  rotor  ie  on  •  vertical  abaft  or  in  a  aero  gravity 
field  the  bearing  load  weald  be  aero  and  tha  eccentricity  ratio 
weald  alee  be  acre.  The  singe  of  tha  shave  carve  M  i  >0  is  the 
film  stiffness  b. 


*.-o  -  o.folaPfi.  ■  »•/*• 
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The  natural  frequency  of  vibration,  atill  considering  tha  ahaft- 
rotor  aystaai  aa  a  simple  rigid  body,  would  than  boi 


*  “  (2°) 

whara  k  -  film  atiffnesa  of  ona  boaring  (lb/in). 

■aa  Figure  (40),  assuning  translation 
only. 

n  ■  total  uaa  of  rotating  system  (shaft) 
(W/g  or  lb.secVln) 

f  “  natural  frequency  (cyclaa/aac) 


In  tha  eaaa  of  this  rotor  tha  waight  ia  0.118  lb. 


0.118 
Tjr  • 


■  •  0.000306  lb.secVln  for  ona  boaring 


For  a  two  boaring  ayataai 

■  ■  0.000306  x  3  lb. sac ^ 

Tha  natural  frequency  than  ia,  fro*  aquation  (20) , 


f  .6x10* 

f  »  964  cyelaa/aoc  or  S7.800  cycles/min. 

Sinca  this  davico  ia  dasignad  to  run  at  24,000  RFN  thara  will 
ba  no  nood  to  paaa  through  this  particular  critical  apaod. 

As  a  further  illustration,  suppose  tha  shaft-rotor  syats*  was 
loaded  by  seaw  force,  so  that  «  was  equal  to  0.5.  A  tangent  to  tha 
curve  in  Figura  (38)  would  yield  a  fila  stiffness  k  ofi 


k  (3.4  -  0,13) lb. 

(0.0002-0.00003)  in. 

k  -  21,800  lbfin. 


OlflMI 


5 


And  the  natural  frequency  corresponding  to  tha  synchronous  resonant 
whirl  or  inversion  point  critical  speed  would  ba  in  this  case  fro* 
aquation  (20), 


f  »  1345  cycles/sec. or  80,600  cycles/fein. 


Although  this  phenomenon  is  rare  in  oil- lubricated  bearings 
because  of  tha  high  dwspinq  in  oil  films,  it  has  been  observed  and 
identified  in  oil- lubricated  shaft  journal  bearings  with  overhung 
compressor  wheels  operating  at  fairly  high  speeds.  Ona  case  involved 
a  shaft  3.3  inches  ia  d lass ter  operating  at  12,000  RFM  with  a  rotor 
weighing  about  100  pounds.  Because  of  the  lew  doping  properties  of 
gas  fllaa  however,  such  instabilities  are  not  suppressed  and  conse¬ 
quently  will  appear  if  tha  dynamic  conditions  are  appropriate. 
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Conical  Synchronous  Whirl  Rotating  Journal 


It  was  assumed  in  the  previous  discussion 
that  the  journal  axis  remains  parallel  to  the  bearing 
axis  as  depicted  on  Figure  (39A).  This  is  termed  the 
cylindrical  mode  or  translatory  mode  of  vibration. 

It  is  also  quite  possible  for  the  shaft  rotor  system 
to  vibrate  in  a  conical  mode  as  illustrated  in  Figure 
( 39B) .  In  this  case  the  journal  axis  generates  a  coni¬ 
cal  surface  with  the  apex  in  the  plane  of  the  center 
of  mass.  The  difference  between  the  two  modes  of  vibra¬ 
tion  is  determined  by  the  distribution  of  the  rotor  mass 
and  for  the  case  of  a  two  bearing  support  system,  the 
distance  between  the  bearings. 


- <£  Bearing 


Fig.  39  -  Indication  of  Synchronous  Whirl  of  Shaft  in  Journal  Bearing. 
(A)  Translatory  or  Cylindrical.  (B)  Conical. 
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In  translatory,  cylindrical,  synchronous  whirl  of 
the  shaft  we  have  the  simple  case  of  an  undeflected  mass 
on  two  springs  with  rigid  bearing  housings.  This  has  al¬ 
ready  been  described.  The  conical  type  of  shaft  synchron¬ 
ous  whirl  is  illustrated  in  Figure  (39)  and  also  schematic¬ 
ally  in  Figure  (40).  As  has  been  explained,  the  springs 
are  not  linear,  as  shown  in  Figure  (38)  but  by  taking  a 
tangent  to  the  load-deflection  curve  they  can  be  assumed  to 
be  linear  at  the  particular  point  of  operation.  For  conical 
whirl,  the  springs  develop  a  restoring  moment,  and  for  two 
springs  this  becomes  kL2/2  ( in. lb. /radian) . 


Fig.  40  •  Schematic  Representation  of  Shaft 
Conical  Synchronous  Whirl. 


Then  the  restoring  spring 
moment  equals  the  inertia 
moment ,  or 


f 


conical 


shaft 


1 

7 


(21) 

where  Its  is  the  trans¬ 
verse  moment  of  inertia 
of  the  shaft  rotor  system 
about  a  diameter  (lb. in. sec2) 


If  in  the  conical  whirling  motion,  the  polar  moment  of 
inertia  of  the  rotating  shaft  system  is  significant  compared 
to  the  transverse  moment  of  inertia,  there  could  be  a  gyro¬ 
scopic  effect  which  would  influence  the  natural  frequency  of 
vibration  and  the  critical  speed  at  which  resonance  takes 
place.  Den  Hartog,  "Mechanical  Vibrations"  4th  Edition, 
McGraw-Hill  Publishing  Company,  shows  that  if  the  shaft  is 
whirling  (or  processing)  in  the  direction  of  spin,  the  gyro¬ 
scopic  inertia  moment  tends  to  make  the  aaiplitude  of  vibra¬ 
tions  smaller,  thus  making  the  effective  spring  rate  higher 
and  in  this  way  raising  the  natural  frequency  of  vibration. 
This  is  the  case  corresponding  to  synchronous  whirl.  Follow¬ 
ing  through  Den  Hartog 'a  dsrivation,  we  have 


^conical  shaft 


(22) 


Where  I_B  is  the  polar  moaient  of  inertia 
of  the  shaft- rotor  aystsm. 
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Synchronous  Trans latorv  and  Conical  Bearing  Whirl 
(a)  Mon- Rotating  Bearing 

Next  we  consider  the  case  of  synchronous  bearing 
whirl.  Such  an  occurrence  is  common  if  the  bearing  is  held 
in  flexible  supports  and  can  vibrate  relative  to  the  spin¬ 
ning  shaft.  We  then  have  the  usual  cases  'f  translatory 
synchronous  whirl, (Fig.  41), 


f 


trans.brg. 


(23) 


and  for  conical  whirl 

^conical  brg. 


(24) 


where  k  is  the  restoring  moment  produced  by  the  gas  film 
in  the  bearing  and  1^  is  the  transverse  moment  of  inertia 
of  the  bearing  about  a  diameter.  An  expression  for  kc  can 
be  obtained  directly  from  the  translatory  k  if  it  is  assumed 
that  the  gas-film  spring  action  in  the  bearing  is  uniform¬ 
ly  distributed,  neglecting  end  effects.  This  assumes  that 
all  sections  of  the  bearing,  including  the  ends,  contri¬ 
bute  equally  to  the  translatory  stiffness.  This  is  not  true 
and  only  in  the  limit  when  1/d  approaches  infinity  would  it 
be  considered  to  be  a  fair  representation  of  the  facts.  How¬ 
ever,  using  this  assumption  and  integrating  across  the  bear¬ 
ing  we  find  that* 


where  1  is  the  length  of  the  bearing.  Now  Equation  (25)  may 
in  general  be  written  as  follows t 


k 


c 


(26) 


where  k*  is  a  numerical  factor  depending  upon  the  1/d  ratio 
of  the  bearing.  Then  Equation  (24)  when  combined  with  (26) 

becomes 
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(27) 


^conical  brg. 


Fig.  41  -  Conical  Whirl  of  Bearing 
With  Respect  To  Shaft. 


1  /  kl* 

2w 


The  number  12  in  Equation 
(25)  and  K*  in  Equation  (26) 
may  be  replaced  by  modified 
values  obtained  from  actual 
experiments  such  as  in  Ref. 
(35).  In  this  paper  exper¬ 
imental  results  are  shown 
for  conical  whirl  of  bear¬ 
ings  (as  contrasted  to  shaft 
whirl)  in  which  this  numeri¬ 
cal  factor  is  evaluated  from 
the  observed  data. 


The  following  values  summarize  the  results t 


TABLE  XVI 


Experimental  Values  of  K*  (Ref. 35) 

tearing  1/d  Ratio  K* 

12  (theoretical  value) 
15 
23 
25 
25 


3 

2 

1.5 

1.0 


Thus,  while  these  values  are  experimental  and  subjected 
to  a  number  of  probable  errors  they  are  nevertheless  more 
accurate  than  the  idealised  value  of  K*  ■  12  computed  for  the 
case  of  1/d  of  •*.  It  is  recommended  that  these  experimental 
values  be  used  in  making  calculations  involving  conical  bear¬ 
ing  whirl  and  in  Equation  (27),  a  number  selected  from  the 
above  table,  for  the  corresponding  1/d  ratio  of  the  bearing 
would  be  used  instead  of  12. 
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(b)  Rotating  Bearing 


Fig.  42  -  Fixed  Shaft  and 
Rotating  Bearing. 


Now  Equation  (27)  for 
conical  whirl  of  the  bearing 
is  applicable  to  the  rocking 
or  oscillating  motion  of  a 
non- rotating  bearing.  If  the 
situation  that  needs  to  be 
evaluated  is  that  of  a  sta¬ 
tionary  shaft  and  a  rotating 
bearing,  gyroscopic  inertia 
effects  of  the  rotating  bear¬ 
ing  would  need  to  be  included 
This  would  correspond  physi¬ 
cally  to  an  electric  motor 
driven  fan,  where  the  Mr«torM 
in  the  classical  sense,  is 
held  stationary,  and  the  hous 
ing  revolves  about  it.  Figure 


(42). 

In  this  case  the  polar  moment  of  inertia  might  be  relative¬ 
ly  large  and  should  be  included.  The  equation  for  natural 
frequency  would  become 


f 


conical  brg. 


i-l/  *  1* 

2*  »k*  (it*-  ipbI 


(28) 


where  I  .  is  the  polar  moment  of  inertia  of  the  bearing  and 
K*  is  a" number  taken  from  Table  XVI  and  which  depends  upon 
the  1/d  ratio  of  the  bearing. 

Let  us  now  summarize  the  discussion  so  far  on  syn¬ 
chronous  resonant  whirl i 


1.  Translatory  synchronous  shaft  whirl  Eq.  (20) 

2.  Conical  synchronous  shaft  whirl 

a.  Without  gyroscopic  effect  Eq.  (21) 

b.  With  gyroscopic  effect  Eq.  (22) 

3.  Translatory  synchronous  bearing  whirl, 

with  rotating  or  non- rotating  bearing  Eq.  (23) 

4.  Conical  synchronous  bearing  whirl 

a.  Non- rotating  bearing  Eq.  (27) 

b.  Rotating  bearing  Eq.  (28) 
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It  is  frequently  possible  to  predict  which  mode  of  syn¬ 
chronous  vibration,  may  be  expected  to  occur  first  as  the 
speed  of  a  unit  is  slowly  raised.  For  the  case  of  shaft 
synchronous  whirl  we  can  take  the  ratio  of  Equation  (20) 
to  Equation  (22)  and 


(29) 


Thus,  if  the  riqht  hand  side  of  the  equation  is  greater 
than  one,  f  ,  conical  shaft  synchronous  whirl  will  occur 
at  a  lower  frequency  than  fts,  translatory  shaft  synchron¬ 
ous  whirl.  If  the  right  hand  side  is  less  than  one  then 
it  follows  that  the  translatory  mode  will  have  the  lower 
frequency. 

In  Ref.  (31)  is  data  for  an  electric  motor  rotor 
held  in  two  rigid  non- rotating  bearings  so  that  Equation 
(29)  would  be  applicable. 


Weight  of  rotor  ■  17.19  lb. 

L  (distance  between  bearings)  -  11.75  inches 
It,  -  1.91  lb. insec2 

I__  ■  0.063  lb.in.sec^ 

P* 
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Substituting  in  Equation  (29) 


mm  ,2 

f£te]  _  4  [1,91-0.063]  _  i  2 

LfcaJ  17.19/386  x  11. 752 


Tha  results  of  tests  made  by  the  authors  of  Ref.  (31), 
although  performed  for  half-frequency  whirl,  did  show 
that  conical  whirl  was  observed  at  a  lower  speed  than 
translatory  whirl,  as  we  would  have  expected  from  the 
ratio  of  ttu2/*cm2  -1.2. 

For  the  case  of  bearing  whirl  (non-rdating)  we  take 
the  ratio  of  Equation  (23)  to  (27),  or 


(30) 


where  values  of  K*  would  be  obtained  from  Table  XVI. 

As  part  of  the  experimental  work  of  Ref.  (35)  it  was 
possible  to  take  a  given  bearing  and  add  weights  to  it  thus 
changing  the  ratio  of  1^  to  m^  in  Equation  (30)  and  in 
turn  changing  the  form  of  the  mode  of  the  lower  frequency 
whirl.  The  observations  were  actually  made  on  half-frequency 
whirl  but  the  general  conclusions  are  valid,  for  synchronous 
whirl,  with  a  non- rotating  bearing. 

Bearing  diameter  -  2  inches 

Bearing  length  -  6  inches 

1/d  ratio  -  3 

K*  from  Table  XVI-  15 
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Therefore,  from  right  hand  side  of  Equation  (30) 


15 

®b  l2 

Mode  of 

Whirl 

Predicted 

Mode  of 
Whirl 
Observed 

1.305 

Conical 

_ 

0.965 

Cylindrical 

Cylindrical 

1.090 

Conical 

Conical 

1.035 

Conical 

Cylindrical 

1.11 

Conical 

Conical 

These  results  show  the  general  agreement  with  the  pre¬ 
dictions  of  Equation  (30). 

For  the  case  of  bearing  whirl,  with  a  stationary 
shaft  but  with  a  rotating  bearing,  we  take  the  ratio  of 
Equation  (23)  to  Equation  (28)  and  havet 


1 -  ]/  *  V? . 

71 II**  litb-  !Pbl 


or 


[L*] 

£cb 


K*  [ltb~  W 


Influence  of  External  Vibration 


(31) 


In  the  previous  section  it  has  been  shown  that 
a  gas-borne  rotating  system  has  at  least  two  resonant  syn¬ 
chronous  speeds  which  correspond  to  natural  frequencies  of 
vibration.  If  such  a  system  is  running  smoothly,  not  at 
one  of  these  "critical  speeds”,  difficulty  may  arise  if  the 
entire  structure  is  shaken  or  vibrated  at  a  frequency  which 
corresponds  to  one  of  these  critical  speeds.  Such  resonant 
external  excitations  are  generally  more  serious  with  gas- 
lubricated  bearings  than  with  liquid  or  oil  lubricated  bear¬ 
ings  because  the  gas  bearing  has  much  less  internal  damping . 
Because  of  this,  it  is  recommended  that  a  bearing  not  be 
subjected  to  a  vibration  excitation  higher  than  75  to  85% 
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of  the  lowest  natural  frequency  of  vibration.  In  the 
case  of  Example  9  for  instance,  a  synchronous  resonant 
shaft  frequency  was  computed  as  964  cycles/sec.  The 
recommendation  then  would  be  to  limit  the  vibrational 
excitation  of  this  unit  to  725  to  825  cycles  per  second 
as  a  maximum. 

If  for  some  reason,  perhaps  test  specifications,  it 
would  be  necessary  to  subject  the  rotor  system  of  Example 
9  to  say  1000  cycles  per  second  it  would  be  desirable  to 
have  the  lowest  synchronous  whirl  frequency  some  15%  to 
25%  above  this  value.  As  shown  in  Equation  (20),  or  sub¬ 
sequent  Equations  (21)  through  (28),  whichever  would  apply, 
this  may  be  accomplished  by  increasing  the  spring  rate  k 
or  kc  and/or  decreasing  the  mass  m  or  the  moment  of  inertia. 

Although  the  damping  in  a  self-acting  bearing  of  this 
type  is  relatively  low,  it  is  possible  to  pass  through  a 
synchronous  resonant  critical  speed,  if  absolutely  neces¬ 
sary,  and  no  damage  may  result  providing  the  amount  of  un¬ 
balance  is  low  and  the  acceleration  rate  is  high.  Under 
no  circumstances  however  should  steady- state  operating 
speed  be  within  approximately  15%  of  one  of  the  natural 
frequencies  of  the  system. 

In  consideration  of  external  vibrations,  resonant 
behavior  of  self-acting  bearings  has  also  been  observed 
at  applied  vibration  frequencies  of  one  half  rotational 
frequency  (Ref.  41).  Of  concern  at  this  condition  is  the 
reduction  in  load  capacity  of  the  bearing.  Observed  in¬ 
creases  in  eccentricity  ratio  from  0.05  at  a  fixed  load 
condition  to  0.6  at  a  sinusoidal  load  condition  and  an 
applied  frequency  of  one  half  rotational  frequency  have 
been  recorded. 

Experiments  conducted  by  Rotron  have  indicated  that 
the  reduction  in  load  capacity  occurs  over  a  rather  narrow 
band  around  the  ratio  of  applied  vibrational  frequency  to 
rotating  bearing  frequency  *  *5.  in  the  case  of  one  parti¬ 
cular  unit,  a  reduction  in  sinusoidal  load  capacity  from 
8  G’s  to  4  G's  was  observed  at  the  applied  frequency  of 
125  cps  with  an  operating  speed  of  250  cps.  This  decrease 
in  load  capacity  occurred  over  the  range  from  100  to  150  cps 
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BEARING  LOAD 
CAPACITY 


applied  frequency.  Figure  (43)  presents  a  dimensionless 
plot  illustrating  the  effect  of  this  behavior. 

When  designing  a  gas  bearing  it  is  quite  evident 
that  the  designer  must  carefully  review  the  vibrational 
requirements  that  the  bearing  must  withstand.  In  addi¬ 
tion  to  avoiding  vibrational  loads  of  the  same  frequency 
as  the  resonant  frequency  of  the  bearing,  load  capacity 
must  be  checked  for  sustaining  vibrational  loads  at  fre¬ 
quencies  within  15*  of  one  half  rotational  frequency. 


APPLIED  VIBRATIONAL  FREQUENCY 
ROTATIONAL  FREQUENCY 


Fig.  43  -  Bearing  Load  Capacity  Under  Sinusoidal  Load  Condition. 
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E.  Half-Frequency  Whirl 

On*  of  the  most  serious  forms  of  instability  en¬ 
countered  in  journal  bearing  operation  is  known  as 
"half- frequency  whirl".  The  phencatenon  is  one  of  self- 
excited  vibration  and  is  characterised  by  having  the 
center  of  the  shaft  orbit  around  the  center  of  the  bear¬ 
ing  at  a  frequency  approximately  equal  to  half  of  the 
spinning  or  rotational  velocity  of  the  shaft.  Under  these 
conditions  Hagg  (24)  has  shown  that  the  capacity  of  the 
bearing  to  support  radial  loads  falls  to  sero.  In  Figure 
(44)  this  half- frequency  whirling  motion  is  depicted  where 
Q  is  approximately  half  of  u.Q  is  the  orbital  frequency 
and  «  is  the  shaft  rotational  frequency. 

The  shaft  system  may  be  stable  as  the  speed  is  in¬ 
creased  until  this  threshold  is  reached.  Crossing  this 
threshold  with  further  increase  in  speed  will  bring  the 
system  into  a  region  of  instability  which  becomes  worm 
violent  as  the  penetration  becomes  deeper  until  inevit¬ 
able  seisure  results.  Unlike  an  ordinary  critical  speed 
the  shaft  cannot  pass  through  this  one  and  attain  a  region 
of  stability  on  the  other  side  at  a  higher  speed  as  with 
synchronous  resonant  whirl.  A  typical  instance  would  be 
to  have  a  bearing  running  smoothly  without  any  difficul¬ 
ties  at  say  40,000  rpm  but  seise  and  fail  completely  at 
43,000  rpm.  Failure  in  most  instances  is  instantaneous 
and  complete  as  the  amplitude  of  vibration  becoaws  equal 
to  the  radial  clearance  in  the  bearing.  It  should  be 
pointed  out  rather  emphatically  that  this  type  of  instabil¬ 
ity  is  not  peculiar  to  gas  bearings  but  actually  is  a  seri¬ 
ous  and  continuing  problem  with  high-speed,  lightly- loaded, 
liquid- lubricated  bearings  as  well. 

The  explanation  of  this  instability  and  possible 
methods  of  predicting  its  occunxnce  are  based  on  an  under¬ 
standing  of  what  steady- state  positions  the  journal  assumes 
in  the  bearing  clearance  with  variations  in  speed  and  load. 
For  a  given  bearing  geometry  these  positions  form  a  path 
which  is  predictable  and  consistent.  This  path  is  the  atti¬ 
tude-eccentricity  locus  as  described  earlier  in  this  manual 
in  Section  n  c.  Referring  to  Figures  (26)  and  (27),  the 
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attitude  angle  is  indicated 
by  the  symbol  ♦  and  for  full 
film,  gas- lubricated  journal 
bearings  may  be  obtained  from 
Figures  (30)  through  (34)  as 
a  function  of  bearing  compres¬ 
sibility  parameter  A. 


Generally  at  light  loads 
and/or  high  speeds,  the  eccen¬ 
tricity  ratio  c  is  small,  ap¬ 
proaching  0,  while  the  atti¬ 
tude  angle  approaches  90°.  For 
heavier  loads  the  eccentricity 
ratio  increases  and  the  center 
of  the  shaft  travels  downward 
along  the  shaft  center  locus 
with  the  attitude  angle  approach¬ 
ing  zero  as  the  eccentricity 
ratio  approaches  one. (Mote  Fig. 26) 

Half-frequency  whirl  is  actually  the  result  of  a  driv¬ 
ing  force  acting  on  the  shaft  causing  it  to  whirl  or  orbit 
in  the  same  direction  as  the  spin  of  the  shaft.  This  driv¬ 
ing  force  is  a  component  of  the  hydrodynamic  load-carrying 
film  force  that  is  developed  in  the  clearance  space  of  the 
bearing.  Thus  in  Figure  (45a)  with  the  load  applied  as 
shown,  the  line  joining  the  center  of  the  bearing  with  the 
center  of  the  shaft  can  be  drawn.  The  components  of  the 
hydrodynamic  film  force  acting  perpendicular  to  this  line 
of  centers  can  be  called  F q.  It  is  the  force  Fg  which  pro¬ 
vides  the  driving  action  and  which  can  result  in  the  unstable 
motion  called  whirl.  In  Figure  (45b),  (c) ,  and  (d)  it  is 
shown  that  the  location  of  a  shaft  center  O'  on  the  attitude 
eccentricity  locus  influences  the  magnitude  of  Fg.  In  Figure 
(45b)  is  depicted  a  position  of  the  shaft  where  Fr  is  small 
and  Fg  ia  rather  large.  Fg  will  therefore  tend  to  cause  a 
high  degree  of  whirl  instability. 

In  Figure  (45c)  it  is  assumed  that  the  shaft  center 
is  now  further  down  the  locus,  perhaps  midway  on  the  arc. 

This  would  result  from  increasing  the  external  load,  for 
example,  or  by  developing  an  internal  load  on  the  journal 
through  the  use  of  grooves  or  slots.  Now  in  Figure  (45c) 


Fig.  44  -  Orbital  Whirling  Velocity  of 
Shaft  Center/1  About  Center  of  Bearing 
Compared  to  Shaft  Spinning  Velocity  uj. 
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the  relationship  between  F  and  f0  shows  that  they  are  both 
of  comparable  magnitude.  If  the  center  of  the  journal  now 
finds  itself  much  further  down  on  the  attitude-eccentricity 
locus  as  shown  in  Figure  (45d) ,  through  modification  of 
either  load  or  geometry,  Fr  becomes  quite  small.  This  means 
that  the  driving  force  Fg  is  becoming  negligible  and  that  for 
these  conditions  the  shaft  would  be  very  stable. 
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Thus  the  angular  position  of  F  and  its  magnitude,  or  con¬ 
versely  the  angular  position  or  F  and  its  magnitude  will 
be  fundamental  parameters  in  establishing  whether  half-fre¬ 
quency  whirl  is  a  possibility. 

Half-Frecruencv  Whirl  Stability  for  Liquid  Lubricated 

Self-Acting  Bearings 

Poritsky  (25),  Boeker  and  Sternlicht  (26),  refer 
to  the  analysis  of  whirl  instability  in  terms  of  these  forces. 
The  equations  of  motion  are  established  and  these  equations 
are  analyzed  for  stability  criteria.  A  relationship  is  de¬ 
veloped  to  indicate  the  threshold  of  whirl  instability  as  a 
function  of  the  mass  of  the  system,  the  spring  constant  of 
the  shaft  itself,  and  the  radial  film  stiffness  of  the  lub¬ 
ricating  film  in  the  bearing  as  determined  by  the  force  Fr 
that  has  just  been  described.  If  Fr  approaches  zero  the 
shaft  will  whirl  at  all  operational  speeds.  An  example  of 
such  a  case  would  be  a  vertical  shaft  supported  by  an  un¬ 
grooved  journal  bearing  where  *  would  equal  to  zero  and  4 
would  be  90°. 

The  expression  derived  is  as  follows,  referring  speci¬ 
fically  to  translatory  half- frequency  whirl  of  the  rotor- 
shaft  systemi 


2_ 

2* 


I 


li  +JL] 


(32) 


where  ffc8  ■  shaft  speed  at  which  half- frequency  transla- 

7"  tory  shaft  whirl  begins  (cycles/sec) . 

m8  *  mass  of  shaft- rotor  system  (lb.sec2/in) 

-  shaft  spring  stiffness  (lb/in) 

K2  *  radial  spring  stiffness  of  the  bearing  film 
for  one  bearing,  obtained  by  plotting  radial 
film  force  Fr  against  film  thickness  and  evalu¬ 
ating  the  slope  of  this  curve  at  the  proper 
eccentricity  ratio  (lb/in) 
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For  our  purposes,  where  the  shaft  system  is  such  that 
operation  is  well  below  the  first  mechanical  critical 
speed  (with  elasticity  of  the  shaft  included) ,  we  can 
properly  assume  that  the  shaft  stiffness  is  much  greater 
than  the  radial  stiffness  of  the  gas  film,  or  K^>>K2  so 
that  in  Equation  (32)  above  1/Ki<<1/K2  and  may  be  dropped. 
Therefore,  Equation  (32)  becomes 


Use  of  Equation  (33)  for  liquid- lubricated  bearings 
has  shown  very  satisfactory  correlation  between  predicted 
whirl  threshold  speeds  and  actual  speeds  at  which  such  in¬ 
stabilities  were  observed  to  begin  (26,27).  More  refined 
analyses  involving  extended  digital  computer  calculations 
are  available  but  the  justification  of  this  extra  involve¬ 
ment  and  complication  has  not  yet  been  established  (28). 

It  is  therefore  recommended  that  Equations  (32)  or  (33)  be 
used  for  design  purposes  with  liquid- lubricated  bearings. 

The  method  provides  a  step-wise  sequence  of  calculations 
(as  shorn  in  Example  10,  Steps  1  through  7)  in  which  a  phy¬ 
sical  understanding  may  be  gained  of  the  phenomenon  known  as 
half- frequency  whirl.  With  this  conceptual  understanding 
there  is  developed  simultaneously  the  knowledge  of  the  in¬ 
fluence  of  the  several  variables  such  as  attitude  angle, 
eccentricity  ratio,  radial  and  tangential  fluid  spring  forces, 
speed,  clearance  ratio  and  length- to- diameter  ratio,  on  the 
threshold  of  ha If- frequency  whirl.  It  follows  that  by  chang¬ 
ing  these  parameters  the  threshold  speed  of  instability  can 
also  be  varied.  In  this  way  the  design  aspects  for  control¬ 
ling  this  phenomenon  are  established. 

This  method  involving  the  use  of  Equations  (32)  or  (33), 
is  also  completely  general  and  applies  to  all  sizes  and 
shapes  of  liquid  lubricated  journal  bearings.  If  stabiliz¬ 
ing  forces  are  introduced,  either  external  to  the  bearing  or 
arising  within  the  bearing  clearance  itself,  the  equations 
can  be  used,  just  as  they  are,  without  modification. 

The  reason  this  method  is  being  given  so  much  emphasis 
even  though  it  was  developed  for  bearings  using  incompressible 
lubricants,  is  that  it  provides  the  first  step  in  developing 
a  design  approach  for  gas- lubricated  bearings.  Actually,  for 
values  of  the  compressibility  number  A  less  than  about  1/3, 
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it  may  be  used  directly  for  gas- lubricated  bearings  pro¬ 
viding  that  the  stability  parameter  u^*  is  large,  that  is 
to  say,  equal  to  or  greater  than  0.5.  This  parameter  is 
described  by  Equation  (34).  For  higher  A  however,  two  al¬ 
ternate  approaches  will  be  suggested. 

Half-Frequency  Translatorv  Whirl  Analysis  for  Small 
Values  of  A 


ihti  calculation  will  involve  the  uao  of 
aquation  (33)  and  with  it  wo  will  attwpt  to  oatabliah  whether 
at  a  given  operating  apeed  the  shaft  will  be  stable  or  not  stable 
to  half-frequency  whirl. 


To  provide  correlation  we  will  use  experimental  data  from 
•ternlicht  and  Winn29.  For  a  bearing  with  the  following  speci¬ 
fications  the  onset  of  half- frequency  shaft  whirl  was  observed 
to  occur  at  11,800  UN. 


Bearing  diameter 
Bearing  length 
Radial  clearance  c 
Viscosity  (air) 


2  inches 
2  inches 
0.001643  Inches 
2.7x10""  lb.sec/ln 
14.7  psla 
11.01  pounds 
11.01  pounds 


Radial  load  on  bearing 
Weight  of  ahaft/bearing 
Shaft  is  horlsontal  and  bearing  loading 
is  due  to  gravity 


We  will  begin  by  asswing  that  the  operating  speed  is 
11,800  UN,  aquation  (33)  will  be  evaluated  and  will  yield  an 
answer  for  the  shaft  speed  *ta/2  In  cycles  per  second. 

If  f.a/2  *  “  11,800  RFM,  this  will  moan  that  the  design 

speed  is  at  the  threshold  of  whirl  instability. 


If  f»a/2  *  60  is  greater  than  11,800  RFM,  let  us  say  for 
example  15,000  RFM,  this  means  that  the  shaft  should  be  free  of 
half- frequency  whirl  at  11,800  RFM  and  thus  is  stable.  It  does 
not  mean  that  15,000  RFM  will  be  a  threshold  speed. 

If  ft|/ 2  x  60  is  less  than  11,800  RFM,  say  5,000  RFM,  the 
indication  is  that  this  design  speed  is  unstable  and  that  the 
shaft  will  be  whirling  at  11,800  RFM  if  it  hasn't  failed  before 
at  a  lower  speed.  It  does  not  naan  that  5,000  RFM  is  a  threshold 

speed. 


Bquation  (33)  is  only  a  stability  indicator.  It  tells  whether 
the  system  is  stable  or  unstable  at  a  pre-selected  speed.  It  does 
not  calculate  a  ■critical*  speed  in  the  usual  sense  unless  the 
assessed  speed  coincides  with  the  calculated  value  from  aquation  (33). 

Actually  it  is  seamwhat  similar  to  a  "go" -"no  go"  gage  estab¬ 
lishing  whether  a  manufactured  part  is  within  acceptable  dimen¬ 
sional  tolerance.  It  either  is  or  it  isn't.  In  the  sans  way, 
aquation  (33)  shews  whether  the  shaft  is  stable  with  respect  to 
half- frequency  whirl  or  it  is  not,  evaluated  of  course  with  refer¬ 
ence  to  an  assumed  operating  speed. 
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»fp  1  -  calculate  * 


Fran  Equation  (1)  A  ■  y^—  (§) 

-  1235  rada/aoc. 

r  ■  1  inch 

a  -  , - 1 - ,)2 

1*7f  1.643x10  J 

A  ■  0.504 

tat  2  -  Kstabllah  Load-Eccentricity  Relationship 

Fro*  Figure  14,  (or  1/d-l,  we  obtain  tha  following 
valuaa  at  *  <  0.S04.  Value  of  d  «  1  x  ^  ■  3x2x14.7  ■  58.8 


c 

W(lb) 

0 

0 

0 

0.1 

0.018 

1.059 

0.2 

0.039 

2.293 

0.4 

0.085 

5.0 

0.6 

0.16 

9.4 

0.8 

0.36 

21.3 

step  3  -  Determine  Attitude  Ancles 

From  Figure  33,  (or  1/d  ■  1,  we  have  the  following  anglea 
at  A  -  0.504 1 


M 


Th«n 


f  ff 

0  82.0 

0.1  81.5 

0.2  81.0 

0.4  76.5 

0.6  65.5 

0.8  41.0 


lt»p  4  -  Determine  the  Component  of  w  acting  alone  the  bine 
of  Centere  of  Bearing  and  Shaft 

Referring  beck  to  Figure  (25),  thie  ie  determined  aa 
W  eoe  4  and  acta  along  the  liieof  centers  00‘.  Thie  is  aleo 
indicated  by  Equation  (16)  and  by  the  force  Fr  in  Figure  (45), 
called  the  radial  force  or  the  force  that  acts  along  tha  line 
of  centers. 

Actually,  steps  (2),  (3)  and  (4)  can  all  be  combined  in 
one  operation  by  making  one  large  table  as  shown  below,  values 
of  hQ  are  obtained  from  Equation  (8). 
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I 

I 

I 

I 


c 

dir. 

h(lb) 

4degreea 

cot  4 

w  coat 
(f,) (lb) 

h  in. 

O 

0 

0 

0 

82.0 

0.1392 

0 

0.001643 

0.1 

0.018 

1 .059 

81.5 

0.1478 

0.1565 

0.001477 

0.2 

0.039 

2.293 

81.0 

0.1564 

0.359 

0.001313 

0.4 

0.08S 

5.0 

?6.S 

0.2334 

1.167 

0.000985 

0.6 

0.16 

9.4 

65.5 

0.4147 

3.90 

0.000657 

0.8 

0.36 

21.3 

41.0 

0.7547 

16.07 

0.000328 

Step  5  -  Date rmlnt  Operating  Eccentricity  Ratio 

The  operating  c  for  thia  bearing  with  N  «  11. OX  pounda 
will  ba  ecnewhara  batwaan  t  ■  0.6  and  t  ■  0.6  judging  by 
tha  abova  tab la.  Fraquantly  it  ia  poaaibla  to  interpolate, 
by  inapaction  but  in  thia  caaa  it  will  ba  nacaaaary  to  plot 
a  curve  of  W  varaua  «  in  order  to  get  a  aufficiantly  accurate 
value. 


i  -  jaaanttMLJSa  to  fgaaUaa  lai 

The  value  of  k2  la  obtained  by  plotting  r_  (which  la 
W  cob  4)  varaua  a  and  hQ  and  datanalnlng  tha  elope  of  thia 
curve  at  tha  oparating  valua  of  (. 
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Slop*  of  eurv*  *t  c  •  0.6S  la  y 


or  Kj  -  *  -  o.oiii?7-4!?00328  “  OiMl? 

Kj  -  24,200  lb/in  at  c  ■  0.63 

8t*p  7  -  calcuUt*  ftt/2  frow  muatlon  (33) 


‘H 

'a 

2 

<sf 


*  Vu?w/««^wh*r*  9 


1 848,000 


9.2*10* 


293  eyclaa/sac. 


386  in/sac* 


or  293  x  60  “  17,600  cyclae  par  ainuta . 


Now  the  measured  value  of  threshold  speed  was  11,900 
cycles  per  minute  so  that  the  predicted  value  based  on 
Equation  (33)  is  higher  than  the  measured  value. 


„  Predicted  threshold  value 
the  ratio  T  •  Measured  threshold  value 


17.600 

lit  Wo 


1.49 


This  is  rather  typical  of  the  cxwparison  between  theoretical 
and  experimental  results  for  shafts  in  360°  gas- lubricated 
journal  bearings.  The  ratio  t" is  largely  a  function  of  A, 
but  not  completely,  as  will  be  shown  below.  However,  as  a 
general  rule,  the  smaller  the  value  of  A  for  a  given  case, 
the  smaller  will  be  the  ratio "cT  This  means  as  we  approach 
the  incompressible  film  condition,  or  as  A  approaches  sero. 
Equation  (33)  becomes  more  accurate. 

For  example,  the  graph  (Figure  46)  shows  the  results  of 
a  number  of  calculations  similar  to  that  of  Example  10.  The 
trend  in  the  relationship  between  the  factor  "<?  and  the  value 
of  A  itf  indicated. 
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0.2  04  0.6  0.6  L0  1.2  I A  16  1.8  2jO 

COMPRESSIBILITY  NUMBER  A 

Fig.  46  -  Ratio  C  Plotted  Against  Compressibility  Number  for  Two  1/d  Ratios. 


Although  Figure  (46)  shows  the  trend  of  actual  ob¬ 
served  threshold  speeds  compared  to  predicted  values  using 
Equation  (33),  there  is  another  parameter  that  is  signifi¬ 
cant  in  establishing  the  stability  characteristics  of  gas- 
lubricated  journal  bearings.  This  is  called  the  stability 
parameter  and  has  been  given  several  forms  depending  upon 
the  investigator  (Refers. 32,  33,  and  34). 

The  form  of  the  stability  parameter  that  appears  most 
useful  from  a  bearing  design  point  of  view  is  «»!*.  This  is 
found  in  Ref. (33)  and  (34) 


(34) 


where  w  is  the  shaft  speed, (radians/sec) 

c  is  the  radial  clearance  in  the  bearing  (in.) 

is  the  mass  per  unit  length  (lb.sec2/in2) 
w  is  the  load  per  length  (lb/in) 

The  values  of  *>i*  can  be  calculated  for  the  data  in 
Table  XVII  using  Equation  (34).  The  results  are  presented 
in  Table  XVIII. 
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Example  lit  Consider  the  bearing  which  is  de¬ 
scribed  as  No.  1  of  Table  XVII,  for  A  «  0.549,  the  mass 
of  the  shaft  is  obtained  from  the  weight  of  the  shaft 
associated  with  one  bearing.  This  weight  was  0.33  pounds. 
The  bearing  was  1/2  inch  long.  The  weight  per  unit  length 
was  thus  0.33  x  2  -  0.66  pounds/inch.  The  mass  per  unit 
length 

2 

0.66  lb. sec 
l“  386  "  in.x  in. 

Mx*  0.00171  lb.sec2/in2 


The  shaft  was  not  horisontal  so  that  the  load  on  the 
bearing  was  less  than  the  weight  of  the  shaft. 

Load  on  bearing  w  -  0.2795  lbs. (given) 

Long  per  unit  length  ■  0.2795  x  2 

«  0.559  lb. /in. 

The  observed  threshold  speed  (or  in  general  the  intended 
operating  speed)  *  584  rads/sec. 


Substituting  in  Equation  (34) 

*  coa  0.000533  x  0.00171 

“1  584  07559 

ej*  -  584  1.63  x  10~6 

-  584  x  1.75  x  10"3 

-  0.745 


Thus  we  have  an  additional  parameter  to  describe 
the  stability  behavior  of  the  bearing. 

Making  similar  calculations,  as  Example  11  above, 
for  the  remainder  of  the  bearings  in  Table  XVII  we  have  the 
following  data,  presented  in  Table  XVIII. 
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TABLE  XVIII 


Stability  Data 


No. 

Source 

Reference 

1 

(30) 

2 

(30) 

h 

3 

(30) 

4 

(30) 

5 

(30) 

h 

6 

(30) 

H 

7 

(31) 

l 

8 

(31) 

l 

9 

(31) 

l 

10 

(31) 

l 

11 

(29) 

i 

12 

(29) 

l 

A 

MQll 

u ,  * 

0.549 

3.02 

X 

0.745 

1.89 

4.94 

0.498 

0.113 

1.44 

1.435 

1.33 

5.07 

0.416 

0.360 

2.35 

0.646 

0.206 

1.94 

0.513 

0.736 

1.965 

1.312 

0.604 

1.97 

1.195 

0.475 

1.55 

1.115 

0.847 

1.985 

1.410 

0.504 

1.49 

2.55 

0.455 

1.39 

2.42 

From  Table  XVI II  it  is  possible  to  cross-plot  and 
average  the  values  of  A  versus  HC"  and  *>]*  for  values  of 
1/d  of  *5  and  1/d  of  1.  Tables  XIX  and  XX  present  the  re¬ 
sults  of  the  cross-plotted  data. 


"C" 


1 

2 

3 

4 

5 


TABLE  XIX 

Cross-Plotted  stability  Data 
vs  "c"  and  for  l/d«*i 

A  A  A 

for  for  for  w^*-l,0 

0  0  0 
0  0.20  0.30 

0  0.44  0.65 

0  0.88  1.30 

0  1.68  2.65 
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TABLE  XX 


Cross-Plotted  Stability  Data 
v|  WC"  and  for  l/d»l 

A  A  A 

"Cw  for  _ for  far  m^-1.5 

10  0  0 

2  0  0.7  0.97 

3  0  1.8  0.20 


The  data  from  Tables  XIX  and  XX  are  plotted  in  Figures 
(47)  and  (48) .  The  original  data  points  from  Table  XVIII 
before  cross-plotting  are  also  shown  to  provide  an  indica¬ 
tion  of  the  averaging  that  was  necessary. 


AWSld-ll1  Aa  u  «fpl.  of  tha  u>  of  thi.  nathod, 
consider  a  cam  takaa  tnm  Reference  (JO).  This  has  pre¬ 
viously  baan  considered  la  Rxanple  11. 

•oaring  laagtl  •  1/1  inch 
•oaring  dims  tor  ■  1  inch 
Radial  eloaraaca  ■  0.000533  lnchos 
Ha. a  Of  rotor  par  bsarlag  ■  0.333/306  -  0.000863 
Ib.aacVin 

Pores  on  bsarlag  -  0.279  iba. 

Observed  whirl  thraahold  apaad  »  91  cyelaa/aac. 
Lubricant  la  air  with  *  -  1.6s  x  10'*  rayna 

Tha  ealculatleos  wars  nada  following  tha  pattarn  usad 
In  Cxanpla  10,  with  tha  following  rosultai 

*  •  0.549 

3^--  0.0300 


<  undor  atatlc  conditions  at  93  cyclas/eec.a  0.4S 
•j-  671  lbs/ln 

f.  ■  fron  0q.  (33)  •  181  cyclas/sac,  pradlctad 
2  thraahold  without  nodlfication. 

New  aodlfy  this  valua  of  pradlctad  threshold  In  accordance 
with  Nquatlon  (34)  and  Plguro  (47)  as  was  dona  In  txanpla  11. 

*1*  '  “^T  <M> 

whara  •  naas/unit  length  •  -  0.001714  lba.sacVln2 


and  •  •  93  <  X  ■  S84  radians/sac 

go  SMi™mzssmi 

1  0.54* 

•  0.?48 
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Fig.  47  -  C  Plotted  AgaimtAond  u>,*  for  l/d=J . 


Fig.  48  -  "C’  Plotted  AgainstAandu>|*  for  l/d=l 
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■oh  *o«*c  Flour*  (47)  tor  1/4  ■  H.  with  »^*  •  0.748 
■ad  *  ■  0.54*  lad  tho  intersection  Hill  occur  vary 
oloooly  to  tho  lia*  uhoro  e  »  3. 

Thor* tor*  dlvld*  281  eycl**/**e  (n*  above,  by  3  tad 
th*  rooult  *111  b* 


-  93.6  cycl**/**c 

Thl*  Hould  lndlcit*  that  th*  h*lt-tr*qu*ncy  Hhlrl  thr**- 
hold  ho*  vary  cloa*  to  th*  operating  «p— d  and  that  ln- 
at ability  nould  b*  *xp*ct*d. 


Half-Frequency,  Translatory  Whirl  of  a  Shaft  in 

Infinitely  Long.  Cylindrical  Journal  Bearings 

For  this  approach  we  utilize  Reference  (34) ,  a 
doctoral  dissertation  by  Professor  V.  castelli  at  Columbia 
Undersity.  This  manuscript  provides  an  analytical  evalua¬ 
tion  of  translatory  half- frequency  whirl  of  shafts  in  cy¬ 
lindrical,  360°  journal  bearings  of  infinite  length.  Al¬ 
though  it  applies  particularly  to  bearings  of  infinite 
length,  it  appears  to  have  some  accuracy  in  predicting  the 
stability  threshold  for  shafts  in  very  long  bearings,  such 
as  for  example  where  l/d  values  may  be  2.5,  3  or  larger. 

It  also  appears  to  be  conservative  for  all  360°  bear¬ 
ings  of  finite  length  providing  that  the  eccentricity  ratio 
for  static  loads,  cc,  is  calculated  for  an  equivalent  bear¬ 
ing  of  infinite  length.  This  can  be  done  through  the  use  of 
Figure  (11)  with  a  unit  load  on  the  bearing, lbs/ in  ,  corre¬ 
sponding  to  that  of  the  finite  length  bearing. 

Thus,  if  such  a  calculation  for  a  bearing  of  finite 
length,  let  us  say  1/d  *  2,  using  the  above  approach  based 
on  infinite  length,  showed  that  for  the  selected  operating 
conditions  the  shaft  was  on  the  threshold  of  instability, 
the  shaft  would  in  all  probability  actually  be  stable  since 
it  has  finite  length  and  not  infinite  length. 

The  finite  length  bearing  has  side  flow  and  the  actual 
value  of  c  would  be  greater  than  *c  as  calculated  for  the 
equivalent  infinite  length  bearing.  There  would  be  a  result¬ 
ant  improvement  in  the  attitude  angle  for  the  finite  length 
bearing  over  the  infinite  length  bearing  and  in  addition 
there  would  probably  be  more  damping  than  for  the  infinite 
case. 


91 


All  o£  these  factors  would  tsnd  to  increass  the 
stability  and  raise  the  threshold  speed  over  that  cal¬ 
culated  for  the  equivalent  bearing  of  infinite  length. 

The  procedure  is  to  make  a  plot  of  A  versus  for 
several  values  of  cQ.  These  curves  represent  the  thres¬ 
hold  of  instability.  is  computed  on  a  static  load  basis. 

A  tabulation  of  the  stability  data  from  Ref.  34  is 
shown  in  Table  XXI. 

A  sample  curve  for  cQ  ■  0.1  is  shown  in  Figure  (49). 
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COMPRESSIBILITY  PARAMETER 


Rg.  49  -  Half  Frequency  Whirl  Stability  Plot  for  ^*0.1 
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TABLE  XXI 


Stability  Data  From  Ref.  34 
for  Half-Frequency.  Translatory 
Shaft  Whirl,  in  360°  cylindrical 


Journal 

Bearings  of 

Infinite  Length 

“l* 

A 

0.1 

0.838 

8.016 

0.1 

0.826 

4.004 

0.1 

0.819 

1.996 

0.1 

0.750 

1.0 

0.2 

0.98 

0.9622 

0.2 

1.09 

1.979 

0.2 

1.06 

4.0126 

0.2 

1.10 

8.058 

0.4 

1.40 

0.8209 

0.4 

1.62 

1.859 

0.4 

1.52 

3.9754 

0.4 

1.52 

8.110 

0.6 

2.14 

0.5418 

0.6 

2.07 

1.460 

0.6 

1.94 

3.643 

0.6 

1.91 

7.781 

0.8 

4.27 

0.2052 

0.8 

4.27 

0.6348 

0.8 

2.68 

2.379 

0.8 

2.32 

6.256 

0.9 

10.05 

0.0722 

0.9 

40.1 

0.2221 

0.9 

13.9 

1.017 

0.9 

3.11 

4.246 
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KARIN*  COMPRESSIBILITY  PARAMETER  A 
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STABILITY  PARAMETER,  U>(+ 


Rg.  50  -  Plot  of  Half-Frequency  Translatory  Whirl  Threshold  for  Infinite  Length 

360*  Journal  Bearing. 
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To  illustrate  the  uaa  of  thia  curve  lat  ua 
calculate  the  predicted  ha If- frequency  whirl  threahold  apaad 
for  a  point  of  aqariMatal  data  fron  laf .  29.  Tha  data  era 
aa  follows: 


d  ■  2  in. 
1  ■  4  in. 
1/d  a  2 


C  •  0.000771  in.  ,,  01 

wgt/brg  a  li.oi  Ibjferq.  w  -  *X*«  2.252  lb/in 
«  a  2.7  n  10“’  rayna  (lb.oec/iri  ) 
f a  a  14.7  paia 
A  a  0.62 


^  "  3a£x4  "  °*00713 


Observed  threahold  apaad  for  half- frequency  whirl  3.200  DM 
or  a  a  339  radiana/aae. 

Calculated  value  of  «j*  fron  aquation  (34) 

*  .  33.jQ*m71L*  9«.997A? 

“I  338|f  2.752 


ax*  -  335 


vr 


995  a  10' 


<t  •  -  0.473 


Sow  to  calculate  a  value  of  i#  for  tha  equivalent  bearing  of 
infinite  length,  turn  to  Figure  (11). 

value  of  for  thia  bearing  is  •  0.0936 


m  Figure  (11)  with  A  »  0.62  and  jIL  «  0.0936,  a  la  very 
eloee  to  0.1.  *  ° 

new  locate  the  inter eaction  of  a,*  and  A  in  Figure  (49), 

It  is  slightly  to  the  left  of  the  threahold  lino  for  the  boaring 
of  infinite  length  and  line  therefore  in  tha  stable  region.  If 
tha  paint  had  fallen  to  the  right  of  thia  line  the  indication 
weald  have  been  for  unstable  operation.  Thus  the  lino  represents 
the  threahold  between  stable  ad  unstable  shaft  rotation. 
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■otic*,  hoMvtr,  that  tha  point  (tea  fcaapla  11,  plottad 
in  Figure  (49),  ia  vary  eloaa  to  tho  threshold  lino  shoving 
that  tho  operating  condition  ia  sassvhat  marginal  aa  Car  aa 
stability  ia  concomad.  Actually,  aa  notod  previously,  this 
data  point  vaa  idantifiod  with  a  measured  spa ad  of  throahold 
inatability. 

Inatoad  of  having  a  separata  stability  plot  for  each  value 
of  *0  ainilar  to  Piguro  (49),  it  ia  convenient  to  eeafeine  all 
data  pointa  from  Table  XXI  on  a  single  graph  as  shown  on  Piguro  (90). 

It  ia  understood  that  only  one  of  those  curves  will  bo  per¬ 
tinent  for  any  particular  stability  evaluation  and  tho  aaaas  con¬ 
clusions  will  bo  valid  aa  outlined  in  laaaglo  lit  namely,  if  tho 
calculated  point  fella  to  tha  loft  of  tha  throahold  lino  tho 
operating  condition  ahould  be  stable,  and  to  tho  contrary,  if  tho 
calculated  point  falla  to  tha  right  of  tho  throahold  lino  operation 
should  be  unstable. 


Summary  of  Discussion  on  Half -Frequency,  Trm  slatory 
Ihaft  Whirl 


The  physical  explanation  for  this  instability  phenome¬ 
non  begins  under  Section  E  and  describes  the  action  of  the 
fluid  film  forces  tending  to  drive  the  shaft  in  an  orbital 
fashion  about  the  center  of  the  bearing. 

A  method  of  analysis  and  prediction  for  rigid-shaft, 
translatory  half-frequency  whirl  is  described  and  presented 
as  Equation  (33).  This  equation  applies  to  liquid  lubricants 
but  has  limited  use  for  gaseous  lubricants  providing  A  is 
small,  perhaps  less  than  1/3,  and  if is  relatively  large, 
perhaps  greater  than  1/2.  The  use  of  Equation  (33)  with  a 
gas  bearing  is  dsmonstrated  in  detail  by  Example  10.  Figure 
(46)  indicates  the  approximate  error  that  may  be  incurred  by 
this  method  in  predicting  threshold  speeds  if  A  grows  in  mag¬ 
nitude. 

In  Equation  (34)  the  stability  parameter  e«*  is  intro¬ 
duced  and  then  the  above  method  of  calculation  for  finite 
length  journal  bearings  is  modified  through  the  use  of  this 
stability  parameter,  with  Example  11  and  Table  XVI I I  showing 
the  end  results.  These  are  shown  graphically  in  Figures  (47) 
and  (48)  for  bearings  of  1/d  ■  and  1/d  ■  1  respectively. 

Use  of  the  graphs  and  the  correction  factor  "C"  for  improved 
accuracy  is  shown  in  Example  12. 
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Finally,  a  method  for  the  prediction  of  half¬ 
frequency,  translatory  whirl  of  a  shaft  in  infinitely 
long,  cylindrical  journal  bearings  for  any  value  of  A  and 
u1*  is  presented  in  Figure  (50).  This  method  can  be 
applied  to  bearings  with  1/d  ratios  of  perhaps  2.5  to  3.0 
and  larger  and  may  also  be  used,  with  some  discretion,  for 
shorter  bearings.  Example  (13)  and  Figure  (49)  demonstrate 
the  procedure. 


Conical  Half-Frequency  Shaft  Whirl  Rotating 
Journal  (Rigid  Body  Condition) 

Conical,  half- frequency  whirl  is  somewhat  simi¬ 
lar  to  the  conical  synchronous  resonant  whirl  described 
earlier  by  Equations  (20)  through  (31)  except  of  course 
that  now  the  orbital  frequency  is  about  half  of  the  spin¬ 
ning  frequency  and  also  the  spring  rate  instead  of  being  k 
as  defined  in  Figure  (38)  is  now  K2  as  defined  in  Eq.  (32). 

Thus,  like  Equation  (21),  we  will  have  for  conical 
half- frequency  shaft  whirl 


L  is  the  distance  between  bearings  as  shown 
in  Figure  (40) 


It8  is  the  transverse  moment  of  inertia  of 
the  shaft  rotor  system  about  a  diameter 
(lb.in.sec^) 

With  conical  motion  of  the  shaft,  there  may  be  gyro¬ 
scopic  moment  contributions  to  increase  the  apparent  stiff¬ 
ness  of  the  shaft,  especially  if  the  mass  is  concentrated 
in  the  shape  of  discs  or  rotors.  Following  the  direction 
of  Den  Hartog  in  "Mechanical  Vibrations",  4th  Edition, 
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McGraw-Hill  Book  Company,  we  arrive  at  Equation  (36)  for 
conical  half- frequency  shaft  whirl. 


l2 

2  *  V2  [Xt8-  2  Ip8] 


(36) 


where  I  is  the  polar  moment  of  inertia  of  the 
P  shaft-rotor  system. 

Conical  Half-Frequency  Bearing  Whirl  (Rotating 
and  Non-Rotating  Bearings) 

Next  we  consider  the  case  of  half- frequency 
bearing  whirl.  This  can  occur  (as  with  synchronous  whirl) 
if  the  bearing  is  held  in  flexible  supports  and  moves  re¬ 
lative  to  the  spinning  shaft.  With  translatory  half-fre¬ 
quency  whirl  of  the  bearing  we  have  Equation  (37) 


(37) 


For  conical  whirl  of  a  non-rotating  bearing  we  have 
Equation  (38) 


where  K2c 


(38) 

(39) 


where  K*  is  a  numerical  factor  determined  from  Table  XVI. 

1  is  the  length  of  the  bearing  as  shown  in  Fig.  (41) 
is  the  transverse  moment  of  inertia  of  the  bearing 
about  a  diameter. 


Therefore,  Equation  (38)  becomes 
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(40) 


a  i IhZ. 

*  V  K*  Itb 


And  finally,  for  the  case  of  a  stationary  shaft  and 
a  rotating  bearing  as  shown  in  Figure  (43) 


f cb  .  1 ■  /  K2  l2 

2  'VK MItb-  2  Ipbl 


(41) 


Let  us  now  summarize  the  various  equations  that 
apply  to  half- frequency  whirl: 

1.  Translatory,  ha If- frequency  shaft  whirl  small 
values  of  A  -  Equation  (33) 

2.  Correction  Factors,  Figures  (47)  and  (48) 

3.  Translatory,  half-frequency  shaft  whirl, 
infinitely  long  bearings  -  Figure  (50). 

4.  Conical,  half-frequency  shaft  whirl 

(a)  Without  gyroscopic  action  -  Eq.  (35) 

(b)  With  gyroscopic  action  -  Eq.  (36) 

5.  Translatory,  half-frequency  bearing  whirl,  with 
rotating  or  non- rotating  bearing  -  Eq.  (37) 

6.  Conical,  half-frequency  bearing  whirl 

(a)  Non- rotating  bearing  -  Eq.  (40) 

(b)  Rotating  bearing  -  Eq.  (41) 

As  with  synchronous  whirl  it  is  frequenly  possible  to 
predict  which  mode  of  half- frequency  vibration  may  be  ex¬ 
pected  to  occur  first  as  the  speed  of  a  unit  is  slowly 
raised.  Far  the  case  of  shaft  half- frequency  whirl  we  may 
take  the  ratio  of  Equation  (33)  to  Equation  (36) 
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Assuming  that  the  value  of  K,  at  translatory  and  conical 
operating  speeds  is  the  same* 


r 

f 


(42) 


Thus,  if  the  right  hand  side  of  Equation  (42)  is  greater 
than  one,  then  f c  /j *  the  conical  shaft  half- frequency 
whirl,  will  occur  fct  a  lower  speed  than  ft  /2,  the  transla¬ 
tory  shaft  half- frequency  whirl.  If  the  right  hand  side  is 
less  than  one  it  follows  that  the  translatory  mode  will  have 
the  lower  frequency. 

Repeating  the  calculations  shown  in  the  section  on 
synchronous  whirl  for  an  electric  motor  rotor  held  in  two 
rigid  non-rotating  bearings  (Ref.  31),  we  have  the  following! 

Weight  of  rotor  •  17.19  lb. 

L  (distance  between  bearings)  -  11.75  inches 
It  «  1.91  lb.  in. sec, 
lp8  "  0.063  lb. in. sec2 

Substituting  in  Equation  (42) , 
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The  results  of  tests  made  by  the  authors  of  Ref.  31 
showed  that  conical  half- frequency  whirl  was  observed  at 
a  lower  speed  them  translatory  whirl  as  would  have  been 
expected  with  a  ratio  of [f ts/2^fcs/2l 2  "  1.16. 

For  the  case  of  a  non- rotating  bearing,  we  take  the 
ratio  of  Equation  (37)  to  Equation  (40). 


and 


(43) 


This  is  the  same  as  Equation  (30)  for  synchronous  whirl 
and  the  experimental  data  taken  from  Ref.  35  establish  the 
validity  of  Equation  (43) . 

For  the  case  of  bearing  half- frequency  whirl  with  a 
stationary  shaft  but  with  a  rotating  bearing  we  take  the 
ratio  of  Equations  (37)  to  (41) . 


yielding 


tb 
_ 2 

fcb 

“2 


k*  [la,-  2  ipb) 
- 


(44) 


102 


lxaaa>ls  14i  Experiment*  at  the  Rotron  Manufacturing  Company 
with  a  stationary  shaft  and  rotating  bearing  of  tha  following  di¬ 
mensions  hava  shown  that  tha  nods  of  half- frequency  whirl  obaarvad 
is  conical.  Let  ua  substitute  tha  particular  values  into  Equation 
(44)  to  saa  if  this  would  hava  bean  predicted.  For  a  hollow  cir¬ 
cular  cylinder  wa  hava  tha  following  relationships.  Tha  symbols 
are  shown  in  Figure  (52)  and  tha  achanatic  representation  in  Fig. (51) i 


Fig.  5 1  •  Schematic  Representation  of  Rotating  Bearing  On  Stationary  Shaft. 


y 


Fig.  52  -  Dimensions  of  Hollow  Cylinder  for  Equations.  (45)  and  (46). 


whore  M  is  the  weight  of  the  bearing. 
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(46) 


yy 


*  l-w  - 


W 


(* 


♦  rJ) 


It  -  0.4S6  in. 
r  -  0.225  in. 
N  -  0.342  lbs. 
1R  *  2. SO  in. 

Substituting  in  aquation  (45) 


»tb 


<>•>42  ^.*0*  (0.4562  ♦  0.2252  j 


-  0.000519  lb. in. sec2 


Next,  substituting  in  aquation  (46) 


2pb  "  [0.4562  ♦  0.22S2] 

I  .  -  0.0001144  lb.in.ssc2 
P* 

Ths  actual  1/d  ratio  for  this  bearing  is  3.33.  Prom  Tabls  XVI 
K*  is  approx iMta.'y  15.  Substituting  into  aquation  (44)  m  have 


‘4 _ 


ISxfo. 000519 


X  1 


-  2  x  0.0001144] 
25 


■otai  Ths  actual  length  of  bearing  insido  of  cylinder  is  1.25 
inches  therefore  1  -  1.25  and  the  1/d  ratio  is  3.33. 


2 

13  x  0.0002502  x  316 
0.342  x  1.25* 


-  3.15 


So  that  ths  tranalatory  node  would  occur  at  a  higher  epeed  than 
the  conical  node.TMa  was  observed  experimentally. 

for  this  bearing  them 


1.775 


This  would  neaa,  f  .  ■  1.7S  x  f 


Seat  1st  us  attsnpt  to  predict  the  actual  whirl  threshold 
speed  for  this  bearing,  using  the  nethod  involving  Figure  (50). 
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«x«a pie  Hi  If  the  bearing  of  Example  14  has  a  diametral 
clearance  of  0.00066  inchee  and  the  actual  bearing  bore  ia  0.375 
inches,  uaing  air  aa  the  lubricant,  what  would  be  the  predicted 
threahold  apeed  for  initiation  of  conical,  half- frequency  whirl? 
For  correlation  with  analytical  prediction,  we  have  a  meaaured 
threahold  of  3900  RFM. 


Let  ua  coneider  a  t  eed  of  7200  EPM.  Aaauming  the  viacoaity 
of  air  aa  2.7  x  10" 5  lba.aec/in  ,  then  from  Equation  (1) 


Therefore! 


*  ■  If  <i>2 

“  ■  I?  “  2"d20*  -  755  rada/aec. 

3>2 

A  .  0.268 


W 

dlP. 


0,342 


0.375  X  1.25  X  14.7 


0.0496 


From  Figure  (11),  for  an  equivalent  bearing  of  infinite  length 
with  A  -  0.268  and  W/dlPa  -  0.0496  we  have  cQ  *0.13. 

In  order  to  make  uae  of  the  atabllity  curvea  Figure  (50),  it 
ia  neceaaary  to  determine  U|*  from  Equation  (34)  where i 


*1 


a 


For  the  caae  of  juat  a  a lap la  gravity  load  of  the  rotor  on 
the  bearing,  Equation  (34)  becomeai 


Therefore 


“1*  "  if  <34A> 

«1*  -  75,^0$?  ■  755  x  0.925  x  10’3 


«j*  ■  0.698 


In  Figure  (50)  there  la  no  etability  curve  for  the  calculated 
t0  of  0.13  and  it  la  neceaaary  therefore  to  eetlmate  the  location 
of  thla  line.  Ne  are  interested  ia  determining  whether  the  value 
o(a*  ■  0.698  and  A  -  0.268  on  Figure  (50)  ia  located  to  the  left 
(stable  region)  or  to  the  right  (unatable  region)  of  the  a.  -  0.13 
atabllity  line.  If  we  plot  the  point  for  A  ■  0.268, a,*  »  0.698,  it 
would  appear  to  fall  on  the  eatiawted  0.13  atabllity  tine.  This 
would  indicate  that  the  threahold  of  traaalatory  bearing  whirl 
would  be  very  close  to  7200  EFN.  From  Example  12,  concerning  this 
bearing,  it  was  determined  that  conical  bearing  whirl  would  be  en¬ 
countered  before  translatory  whirl  with  the  following  relationship! 
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1.78 


Tharafora  tha  thraahold  of  half-fraquancy  conical  whirl  would  ba 
pradictad  asi 


£  .  7200 

1.78 


4050  RPK 


Aa  wma  atatad,  tha  aaaaurad  valua  waa  3800  RPM. 
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III.  SELF-ACTING  THRUST  BEARINGS 


Self-action  thrust  bearings  may  be  of  many  types. 
All  of  the  usual  varieties  may  be  employed  with  gaseous 
lubricants.  These  are  the  classical,  pivoted-pad  type, 
the  tapered  land  with  fixed  geometry,  the  Rayleigh  step 
bearing  and  various  forms  of  the  pocket  bearing.  Figure 
(53). 


u 


In  addition,  the  spiral 
or  herringbone  grooved 
pumping  plate  may  be  used, 
which  is  in  essence  a  mod¬ 
ification  of  the  Rayleigh 
step  as  applied  to  circular 
geometry.  See  Figure  (54) 
and  also  Figure  (5),  which 
is  reproduced  for  conven¬ 
ience  as  Figure  (55). 


The  choice  among  these 
bearings  will  depend  upon 
the  relative  ease  of  phy¬ 
sical  construction,  prac¬ 
ticability  of  installation 
in  a  given  machine,  and 
ability  to  carry  the  speci¬ 
fied  loads  within  the  limits 
of  available  space.  The 
tilting-pad  type  of  bearing 
is  self-aligning  but  does 
not  equalize  the  loads  on 
each  shoe  of  the  bearing. 

The  other  types  are  not  self-aligning  and  usually  require  some 
type  of  gimbal,  pivoted  or  flexible  mounting  to  perform  this 
function.  Many  of  the  comments  previously  made  for  self-act¬ 
ing  journal  bearings  apply  equally  as  well  to  thrust  bearings. 
Pressures  are  developed  in  the  same  fashion  and  a  typical  pat¬ 
tern  for  a  tilting  pad  type  is  shown  in  Figure  (56).  in  a 
qualitative  way  Figure  (56)  shows  the  difference  between  the 
pressure  rise  for  a  compressible  and  an  incompressible  lubri¬ 
cant.  Figure  (57)  shows  quantitative  values  for  such  a  bearing. 


(C) 


»> 


/Z4Z/////Z//ZZ/Z 


'///Zz 

1 _ 


RoyMgS  Step 


Pocket  Typo 


Fig.  53  -  Typical  Forms  of  Self-Acting 
Thrust  Bearing.  (Ref.  21). 
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Fig.  54  -  Grooved  Pumping  Plate  With 
Either  Logarithmic  or  Archimedean  Spiral. 


Fig.  55  -  Grooved  Thrust  Bearings. 


KARINS  LENSTH 


Fig.  56  •  Comparative  Pressure  Distribution 
for  Compressible  and 
Incompressible  Lubrication. 

(From  Gross,  Ref.  36). 


Note  that  for  high  values  of 
A  (see  Eq.2)(or  in  this  case 
as  the  ambient  pressure  falls) 
the  center  of  pressure  moves 
toward  the  trailing  edge  or 
the  back  of  the  pad.  For  each 
such  operating  condition  the 
pivot  must  be  located  at  the 
center  of  pressure.  Where  the 
ambient  pressure  or  speed  are 
likely  to  change  considerably 
( A  varies  over  a  fairly  wide 
range)  this  would  be  a  disadvan¬ 
tage. 

Thrust  shoes  made  of  fixed 
tapered  lands,  Rayleigh  steps 
or  pockets,  would  not  be  sensi¬ 
tive  to  this  situation. 

The  literature  is  quite  com¬ 
plete  on  the  analysis  and  design 
of  the  tilting-pad  thrust  bear¬ 
ings  (Ref.  11).  However,  they 
will  not  be  considered  here  be¬ 
cause  of  the  practical  limita¬ 
tions  imposed  by  the  small  size 
of  the  blowers  toward  which  this 
design  manual  is  pointed.  A  lower 
limit  on  the  size  of  a  tilting- 
pad  thrust  bearing  with  an  equal¬ 
izer  system  would  be  about 
1>S  inches  outside  diameter.  They 
can,  of  course,  be  made  smaller 
but  the  cost  of  manufacture  would 
make  them  impractical  for  use  in 
commercially  competitive  equipstent. 

Instead,  design  procedures  for 
the  spiral  grooved  pumping  plate 
will  be  described  as  shown  in  its 
optimised  form  in  Figure  (55b). 
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1  Pa*  142.8  MIA,  A*  2 

2  Pa*  14.29  MIA,  A-  20 

3Po*  3.35  P9IA.A*  93 

4 Pa*  1.425  PSIA.A*  200 


Fig.  57  -  Gas  Pressure  Distribution  in  Flat  Tilted  Pad  Bearing. 


Thxs  is  baaed  on  the  principle  of  the  Rayleigh  step  bear¬ 
ing,  Figure  (53c),  where  gas  is  dragged  into  a  slot  by  a 
moving  runner.  The  exit  end  of  the  slot  has  a  restrictor 
or  dam  so  that  the  escape  of  gas  from  the  slot  is  retarded. 


PUMPINQ  ACTION 

Fig.  58  -  Schematic  Representation  of 
a  Pumping  Groove. 


The  runner  continues  to  pump 
gas  into  the  entrance  of  the 
slot  due  to  viscous  drag  and 
as  a  consequence  the  pressure 
in  the  slot  builds  up  and  the 
bearing  has  the  capacity  to 
carry  load  with  a  gas  film 
separating  the  surface.  Fig. (58). 
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Reasonably  high  pressures  can  be  generated  and  since 
the  surface  area  on  such  a  pumping  plate-thrust  bearing 
is  large,  the  load-carrying  capacity  is  surprisingly  high. 

A  detailed  hydrodynamic  theory  for  the  spiral  grooved 
thrust  plate  has  been  worked  out  by  Whipple  (37).  This 
theory  includes  the  types  shown  in  Figure  (55a)  and  (b) .  He 
has  shown  that  for  the  herringbone  thrust  plate  the  pressure 
rise  developed  by  the  grooving  is  independent  of  compressi¬ 
bility  effects.  For  the  spiral  grooved  thrust  plate  Whipple 
indicates  that  compressibility  effects  do  occur  but  the  ef¬ 
fects  are  not  significant  until  the  pressure  rise  developed 
by  the  grooves  is  three  or  four  times  the  ambient  pressure. 
For  applications  in  normal  atmospheres  such  limits  are  not 
generally  exceeded  and  as  a  consequence  the  solution  of  the 
equations  based  on  an  incompressible  lubricant  may  be  used. 
This  may  be  measured  in  terms  of  another  form  of  the  compres 
sibility  bearing  parameter  A  .  For  thrust  plates  A  may  be  de 
scribed  as* 

A  =^22  (3) 

h  Pa 

where  J  is  the  radial  dimension  of  the  thrust  plate. 

The  data  show  that  for  values  of  A  up  to  30  very  little 
difference  is  discernable  between  solutions  for  incompress¬ 
ible  and  compressible  lubricants  as  applied  to  the  spiral 
grooved  pumping  plate. 

Or.  Boeker  (38),  reports  on  a  detailed  calculation  for 
a  grooved  thrust  plate  of  this  type  having  an  outside  dia¬ 
meter  of  four  inches  and  an  inside  diameter  of  2.5  inches, 
operating  in  air  at  atmospheric  pressure.  With  a  running 
film  thickness  of  0.0003  inch  and  a  speed  of  10,000  RPM  and 
using  the  optimised  geometry  recommended  by  Whipple,  the 
average  film  pressure  developed  is  about  4.6  psi.  With  an 
effective  plate  area  of  7.65  square  inches,  the  load-carry¬ 
ing  capacity  is  calculated  to  be  about  35  pounds. 

Whitley  and  Williams  (39)  indicate  however,  that  unless 
the  plates  of  the  thrust  bearing  are  especially  rigid,  dis¬ 
tortions  and  deflections  of  the  plates  will  reduce  the  esti¬ 
mated  load-carrying  capacity.  Depending  upon  the  degree  of 
distortion  the  measured  load  has  ranged  from  *3  to  3/4  of  the 
predicted  theoretical  load.  Other  investigators  show  results 
ranging  mostly  from  3/4  up  to  the  full  theoretical  value  of 
load. 
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Spiral  Grooved  Thrust  Bearings 

A  single  spiral  groove,  as  shown  in  Figure  (54),  will 
produce  a  pressure  rise  and  some  load-carrying  capacity. 
However,  if  the  design  is  optimized  for  maximum  load-carry¬ 
ing  ability  a  multiplicity  of  shorter  grooves  will  be  in¬ 
dicated.  This  has  been  done  by  Whipple  (37). 

He  concludes  that  there  should  be  18  grooves  of  width 
a^.  The  land  width  between  grooves  is  a2.  See  Figure  (59). 


Dimension  ct  is  an  ungrooved  annulus  known  as  the 
seal  belt.  This  corresponds  to  the  dam  or  flow  restric¬ 
tor  in  the  Rayleigh  step  bearing.  Figure  (58).  In  this 
case,  shown  in  Figure  (59),  the  bearing  is  considered  an 
inward  pumping  surface  with  the  direction  of  rotation  of 
the  mating  thrust  plate  CCW.  If  the  direction  of  rota¬ 
tion  is  reversed  and  the  seal  belt  located  toward  the  outer 
diameter  of  the  thrust  plate  we  would  obtain  an  outward 
pumping  thrust  surface. 

In  Figure  (59) t 

(6 . )  is  the  depth  of  the  groove 
(h)  is  the  film  thickness  of  the  bearing  mea¬ 
sured  between  the  top  of  the  land  and  the 
thrust  face. 

Whipple  has  established  that  for  optimum  load-carrying 
capacity  the  following  relationships  should  be  satisfied. 


3.05 

(47) 

1.8  a2 

(48) 

0.7 

(49) 

72° 

Slight  deviations  from  the  optimized  parameters  are 
possible  without  too  significant  a  change  in  performance. 
For  greater  variations  actual  evaluations  by  testing  is 
recommended. 


Two  parameters  for  the  analysis  of  these  optimized 
bearings  have  been  developed  by  Whipple.  They  are* 


e 


T 


T 

(b+c)  Pa 


1 

1.145 


X  uQ(b+c) 


P 

a 


(50) 

(51) 
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VALUES  of  9 


where  the  individual  variables  are  described  as  follows: 


T  *  thrust  capacity  per  inch  of  mean  circum¬ 
ference  (lb/in) 

u  *  the  linear  velocity  at  the  mean  radius 
(in/sec) 

Values  of  t  and  6  are  given  in  Table  XXII. 

TABLE  XXII 


Values  of  Spiral-Grooved  Thrust 
Plate  Parameters,  0  and  t  (Ref. 37) 


T 

10 

8 

6 

4 

2 

0 


_G _ 

2.62 

2.25 
1.78 

1.26 
0.70 

0 


VALUES  of  r 

fig.  60  -  Spiral-Grooved  Thrust  Plate  Parameters  from  Table  20.  (37). 
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Example  l6i  To  illustrate  the  design  approach  established 
by  Nhippla  lat  ua  determine  tha  load-carrying  capacity  of  a  par¬ 
ticular  thruat  boaring  geometry  having  an  outaida  diameter  of 
1.365  inchaa  and  an  inaida  diameter  of  0.420  inchaa.  Rotor  apaad 
3450  RPM.  Rotor  weight  0.388  lba.  Magnetic  thruat  load  0.30  lba. 
Referring  to  Figure  (59), 


-  0.4 20f|* 

r 

> 

V 

J 

ct+  b 


1.365-0.420 


-  0.473" 


1.365“ 


From  Equation  (49) , 


b  -  0.7 


<h+c} 


0.7(0.473)  •  0.3311" 


or  ct  -  0.473-0.3311  -  0.142" 

Mean  radio,  r  -  -  2*^2 

m  - n -  +  0.210 


r  -  0.446" 

m 


Mean  Circumference  ■  2*  r^ 

Mean  circumference  -  2*  x  0.446 

Mean  circumference  ■  2.80" 

2  80 

With  18  groovea,  a^  ♦  e^  ■  — Jg- 

al  +  a2  •  0.1555" 


From  Equation  (48)  1.8 
Subatituting  in  above 

*2 

and  a 

1 


•l 


1.8  »2  ♦  *2 

2.8  a. 


0.1555" 

0.1555* 

0.0556’ 

0.0999" 


•  g 

Now  with  air  having  a  viscoaity  of  2.86  x  10  reyne,  we  will 
next  evaluate  the  parameter  t.  The  peripheral  speed  at  the  aMan 
radius  is 


3450 


x  2s  x  0.446  •  161  in/sec. 
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Than  In  Equation  (Si) 


nu. 


■  1.145 


<bn> 


.  1 

1.145 

.  12*28 
it 


3.86xl0~9x  161  g  0.473 
h*  *  14.? 


«  10 


-9 


aaaua 

than 


h  -  0.0001* 

ia.»e  x  10* 9  12.98  x  10'9 


*  “  - XT 

(1x10  ’) 4 


10- • 


1.398 


and  from  rigura  (60)  8  ■  0.454 


For  aatall  valuaa  of  t  and  8  straight  lina  interpolation  My  ba 
Mda  with  good  accuracy  ainca  tha  plot  in  Figure  (60)  ia  essen¬ 
tially  linear  for  this  range  of  variables.  For  example. 


X  0.7 
1.298  ”  2.0 


*  “  fto  *  K29e 


x  -  0.454  ■  8 


Than  f ran  Equation  (50) t 


8  - 


(b~t>F. 


t  -  8  p>*ct)  rm 

r  -  0.454  X  0.473  x  14.7 

T  ■  3.16  lbs/in . thrust  capacity  par 
inch  of  assn  circumference 

Total  thrust  capacity  ■  T  x  2*rB 

-  3.16  x  2.80 

-  8.85  lbs  with  h  -  0.0001* 

If  this  war  a  tha  appliad  thrust  load,  than  tha  groove  depth 
would  bo  data  rained  by  Equation  (47), 


4t«  3.05  h 

in  this  caaa  *t“  3.05  x  0.0001 

4t-  0.0003  inches 


However,  it  ia  desirable  to  develop  a  aore  coop  lata  picture 
of  tha  bearing  performance.  This  can  be  done  by  assuming  different 
values  fcr  h  and  confuting  the  corresponding  values  of  t  ,  8,  T  and 
total  thrust  capacity,  as  shown  in  Table  XXIII. 
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0.0001 

1.298 

0.454 

3.16 

8.85 

0.0002 

0.3245 

0.1137 

0.792 

2.215 

0.0003 

0.1443 

0.0505 

0.353 

0.984 

0.0004 

0.0811 

0.0284 

0.198 

0.554 

0.0005 

0.0519 

0.01815 

0.1265 

0.354 

0.0006 

0.0360 

0.0126 

0.0878 

0.246 

0.0007 

0.0265 

0.00926 

0.0645 

0.1807 

0.0008 

0.0203 

0.0071 

0.0495 

0.1385 

0.0009 

0.01603 

0.00561 

0.0391 

0.1095 

Figure  (61) 

depicts  a 

plot  of 

the  total 

load-carrying 

city  of  auch  •  thrust  bearing.  On  the  mm  plot  is  also  shown  a 
superposition  of  two  of  these  curves,  bade  to  back,  as  would  be 
the  case  with  a  double  acting  thruet  bearing  having  a  total  clear¬ 
ance  (h,  +  hj)  of  0.001  inches.  The  configuration  of  such  a  double 
acting  beerLvi  is  shown  schsMtically  in  Pi  jura  (62) . 

consider  first  that  we  have  only  a  single  thrust  bearing, 
the  left  hand  bearing.  With  a  thrust  load  of  0.3  pounds  the 
operating  f ila  thickness  would  be  about  0.0005  inehee. 

With  a  double  acting  bearing  and  the  sane  antemal  thrust 
load  the  operating  f ila  thickness  would  be  a  little  lees  because 
the  second  thrust  bearing,  while  not  carrying  any  internal  load, 
would  nevertheless  (npnas  an  additional  force  upon  the  thrust 
face  that  is  already  carrying  load.  The  magnitude  of  this  addi¬ 
tional  internal  thrust  load  will  depend  upon  the  total  clearance 
designed  into  the  double  acting  thrust  bearing.  This  is  the  clear¬ 
ance  deeignated  as  h,  ♦  hj  inehee.  Figure  (62) .  If  a  bearing  is 
built  with  a  total  clearance  of  0.001  inches  ss  depicted  in  Figure 
(61),  then  the  net  load-carrying  load-capacity  of  the  left  hand 
bearing  would  be  given  by  Table  XXIV. 

Table  XXIV 


h 

inches 

T  total 
baft  brg. 
lb. 

T  total 
Right  brg. 
lb. 

T 

Wet  lbs. 

0.0001 

8.85 

0.1095 

8.74 

0.0002 

2.215 

0.1385 

2.076 

0.0003 

0.984 

0.1807 

0.8033 

0.0004 

0.554 

0.246 

0.308 

0.0005 

0.354 

0.354 

0 

0.0006 

0.246 

0.554 

-0.308 

0.0007 

0.1807 

0.984 

-0.8033 

0.0008 

0.1385 

2.215 

-2.076 

0.0009 

0.1095 

8.85 

-8.74 

These  values  are  shown  in  Figure  (63) . 

It  should  be  noted  that  whan  the  runner  is  centered  in  the 
double  acting  thrust  bearing,  correspond ing  to  the  esse  of  h^»  h« 
in  Figure  (62) ,  the  load  capacity  of  tho  thrust  bearing  is  sero 
as  illustrated  on  Figaro  (63)  for  a  file  thickness  of  0.0005  inches. 

If  there  is  a  thrust  force  of  0.3  pounds  acting  to  the  left, 
in  Figure  (63) ,  the  f ila  thickness  for  the  double  acting  thrust  bear¬ 
ing  is  0.0004  inches.  For  the  single  acting  thrust  bearing  with  the 
mow  magnitude  of  load  the  filn  thickneea  would  be  0.0005  inches,  as 
was  pointed  out  previously  in  conjunction  with  Figure  (61). 
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THRUST  LOAD  LBS 
LEFT  HAND  SEARINS 


I - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 

.001  JOOOO  .0000  .0004  ,ooo>  0 

FILM  THICKNESS,  RISHT  HAND  SEARINS  (INCHES) 

Fig.  61  -  Thrust  Bearing  Capacity  of  a  Single  Bearing  and  Also  of  a  Double  Acting 
Bearing  with  a  Total  Clearance  of  0.001  Inches 


Thruot 


Fig.  62  -  Schematic  Diagram  of  Double  Acting  Thrust  Bearing. 
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Fig.  63  •  Thrust  Load  Characteristics  of  Double  Acting  Bearing  With  Total 
Clearance  of  0.001  Inches. 
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NET  THRUST  LOAD  LBS 


Film  Stiffness  in  Spiral  Grooved  Thrust  Bearings 


Where  the  dynamics  of  the  thrust  bearing-rotor 
assembly  must  be  considered,  it  is  necessary  to  evaluate 
the  stiffness  of  the  gas  film  in  the  bearing.  This  stiff¬ 
ness  is  defined  as  the  rate  of  change  of  thrust  capacity 
with  respect  to  change  in  film  thickness  or  dT/dh,  where  T 
is  the  thrust  load  in  pounds  and  h  is  the  film  thickness  in 
inches.  The  units  of  dT/dh  are  lbs/in.  so  that  we  are  deal¬ 
ing  with  the  property  of  the  gas  film  that  behaves  like  a 
spring.  This  characteristic  was  evaluated  earlier  for  the 
journal  bearing  when  synchronous  whirl  in  its  various  forms 
was  investigated. 

Referring  back  to  Example  14  and  Figure  (61) ,  with  a 
thrust  load  of  0.3  pounds  the  expected  film  thickness  will 
be  0.0005  inches.  If  we  draw  a  tangent  to  the  load-film 
thickness  curve  at  h  *  0.0005  inches,  the  slope  of  the  curve 
will  be  a  quasi-static  version  of  the  film  stiffness  which 
should  be  satisfactory  for  the  determination  of  the  dynamic 
characteristics  of  the  system  in  the  axial  direction.  For 
example,  the  slope  of  the  tangent  in  Figure  (61)  at  h  ■ 
0.0005  in.  is 

1.15-0.354  _  0.796 
0.0005  0.0005 

therefore  0.796  .. 

*  "  0.0005  "  lba/rn. 


The  natural  frequency  of  axial  vibration  can  then  be 
estimated,  knowing  that  weight  of  the  rotor  is  0.388  lbs., 
by  using  Equation  (20)  after  substituting  the  appropriate 
terms. 


Thus, 


f  1_  J 1592  x  386 

2w  Y  0.388 


f  ■  200  cycles/sec.  (axial  natural 

frequency  for  single-acting  bearing) 


Now  consider  the  double-acting  thrust  bearing  for  the 
same  loading  and  mass  condition.  A  tangent  will  be  drawn 
to  the  curve  at  a  film  thickness  of  0.0004  inches.  The  slope 
of  this  tangent  becomes 


1.17 

0.0005 


3400  lbs. /in. 


1)9 


I 


It  becomes  evident  that  the  effect  of  the  additional 
thrust  bearing  is  to  apply  what  might  be  termed  a  "preload" 
to  the  unit  and  in  this  way  make  the  film  stiffer  and  thus 
raise  the  natural  frequency  of  vibration  of  the  rotor  sys¬ 
tem. 


It  follows  then  that: 

1_  13400  x  38& 

"  2»  y  0.388 

f  *  292  cycles/sec. (axial  natural  fre¬ 
quency  for  double-acting  thrust  bearing.) 

Using  a  total  clearance  less  than  0.001  inch  would  in¬ 
crease  the  magnitude  of  the  "preload"  and  thus  raise  the 
natural  frequency  of  vibration. 

The  most  complete  collection  of  references  on  all 
phases  of  gas- lubricated  analysis  and  design  is  Ref.  (40), 
which  contains  464  abstracts  of  papers  in  this  field. 
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IV.  BEARING  LOADS 


The  determination  of  bearing  loads  is  an  important 
aspect  of  gas  bearing  design.  The  rotating  mass  load  and 
mass  unbalance  load  are  readily  attainable  by  calculation. 
However,  it  is  also  necessary  to  consider  radial  loads  pro 
duced  by  the  motor  itself.  Of  primary  concern  is  the  mag¬ 
netic  force  produced  by  an  eccentric  air  gap  between  the 
motor  stator  and  the  motor  rotor. 

A.  Magnetic  Radial  Bearing  Loads 


> 


In  either  of  the  three  cases  depicted  in  Figure  (64) 
the  eccentricity  produces  an  unequal  air  gap  between  the 
motor  stator  and  rotor  which  results  in  a  magnetic  radial 
force.  This  force  is  due  to  an  unbalance  in  the  radial 
direction  of  the  linking  of  flux  lines  in  the  motor  gap 
region.  In  the  case  of  a  two  pole  motor,  the  flux  lines 
link  the  stator  north  and  south  poles  through  the  rotor  and 
the  stator  as  shown  on  Figure  (65) .  The  magnetic  field  ro¬ 
tates  at  a  frequency  equal  to  line  frequency.  At  the  in¬ 
stant  of  time  when  the  stator  poles  are  in  position  (65a) 
there  is  equal  linking  of  the  flux  lines  through  the  rotor 
and  no  radial  forces  exist.  As  the  field  rotates  to  posi¬ 
tion  (65b) ,  due  to  the  eccentric  rotor  condition,  there  is 
an  unbalance  of  flux  lines  through  the  rotor  producing  un¬ 
balanced  radial  magnetic  force  in  the  direction  indicated  by 
the  vector  F.  The  force  drops  to  zero  with  an  additional  90° 
rotation  of  the  field  (65c)  and  peaks  to  a  maximum  again  at 
270°  total  rotation  (65d) .  it  can  be  seen,  therefore,  that 
the  magnetic  radial  force  is  a  sinusoidal  directional  force 
acting  at  the  minimum  motor  air  gap  region  at  a  frequency 
twice  line  frequency.  The  magnitude  of  the  force  is  a  func¬ 
tion  of  the  eccentricity  and  can  be  described  analytically  as: 


F 


17.42  x  (jqqoo^  x  D  x  l^x  o  (52) 


Fig.  65  •  Flux  Line*  Unking  Motor  Rotor  and  Stator  for  Angular  Position*  of 

Stator  Magnetic  Field. 
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where 


F  =  Force  (peak  value)-lbs. 

B  =  Flux  Density  -lines  per  sq. in. 
D  “  Rotor  Diasieter  -  inch. 


*m 


o 

Y 


A 


Rotor  Stack  Length  -  inch 

f-3  [  ■  2~^-- 1-  1  (53) 

2*  2  1-r2 

eccentricity.  _  ..  . 

- -  (eccentricity  =  2  x  distance 

between  rotor  and  stator 

centers)  (54) 

Average  motor  radial  air  gap  -  inch. 


The  force  described  by  this  forsmla  is  a  peak  value  vary¬ 
ing  from  zero  to  a  saxiww  sinusoidally.  The  average  force 
acting  on  the  rotor  is  one  half  the  peak  value  (F/2) . 

If  the  case  of  Figure  (64b)  is  considered,  i.e.  the 
rotor  rotating  about  a  center  other  than  it's  geometrical 
center,  then  the  location  of  ■inisua  air  gap  rotates  about 
the  stator  at  rotor  rotational  frequency.  In  an  induction 
motor  the  rotor  rotates  at  a  frequency  less  than  the  stator 
field  rotational  frequency.  The  difference  between  stator 
and  rotor  frequencies  is  the  slip  frequency  of  the  notor. 

As  seen  on  Figure  (65) ,  the  radial  magnetic  force  occurs 
when  the  stator  poles  are  at  right  angles  to  the  point  of 
mininun  motor  air  gap.  This  point  of  minimum  motor  air  gap 
is  rotating  at  slip  frequency  in  relation  to  the  stator  field 
rotation,  and  therefore  the  radial  magnetic  force  rotates  at 
slip  frequency. 


The  most  common  type  of  eccentricity  encountered  in 
electrical  motors  is  a  combination  of  both  displaced  rotor 
and  stator  centerlines  and  rotation  of  the  rotor  about  a 
center  other  than  its  geometrical  center.  The  resultant  force 
varies  in  magnitude  at  slip  frequencies  from  a  maximum  when 
both  eccentricity  effects  are  additive  to  a  minimum  when 
both  effects  are  subtractive. 
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The  significance  of  these  forces  can  be  shown  in  the 
following  example  for  a  2  pole.  60  cycle  induction  motors 

Motor  Length  (L)  *  1.250  inch 
Motor  Dianeter  (D)  «  0.923  inch 
Flux  Density  (B)  *  25,000  lines  per  ■Biwthl 
Average  Motor  Air  Gap  (1)  *  .0075  inch 
Eccentricity  due  to  misalignment  of  rotor  and 
stator  centerlines  *  .001  TIB  (inch) 

Runout  of  Rotor  O.D.  to  Bore  *  .0005  TIB  (inch) 
Total  Eccentricity  *  0.001  *  0.005  *  0.0015 


lZi 


(tt»  r  -  17.«  f||i*ir*-  •»  V** 

w  -  H®5  -  ftrSkr  -  •-»  ■*- 

-*»• 

r9mM  -  i7.e  *  x-no  *  •>«» 

r  **«•?»  *  1  -1'* 

T _  -  0.799  rmmtm  (fHkl 


r. 

»^.r 


-  0.797 

-  0.2*1 
-  0.131 


■who  tore*  tfcoicfoc*  to  iv  tW 
9*0  oceerie*  M  120  cycle* 
fore*  lml  of  0.797 
•lip  IraoMKy  of  to* 


124 


gearing^  Wtanfwfe  toads 


Another  condition  that  can  be  oncountarad  in 
electrical  motors  is  tha  non-parallelism  of  tha  stator 

rotor  center lines .  As  shown  on  Figure  (66)  this  con¬ 
dition  creates  am  eccentric  air  gap  at  one  end  of  tha 
rotor  180  degrees  out  of  phase  with  tha  eccentricity  at 
tha  other  end  of  the  rotor,,  thus  producing  a  moment  on  the 
rotor.  As  can  be  seen  from  Equation  (52)),,  the  magnetic 

force  due  to  motor  air  gap  eccen¬ 
tricity  varies  directly  with  both 
tha  rotor  length  and  tha  parameter  a 
For  low  eccentricities,,  at  least 
below  .004  inch  TTH„<r  can  be  as¬ 
sumed  to  vary  directly  with  the 
eccentricity .  Table  1  presents 
tha  magnetic  force  versus  eccen¬ 
tricity  for  Example  15  and  indi¬ 
cates  very  slight  deviation  from 
a  straight  line  relationship  be¬ 
tween  force  aid  eccentricity  to 
eccentricities  of  .004  inches. 


HfliT.I  I 

Eccentricity  of 


motor  Air  (Sag 

T.I -A-Cim* 

Magnetic  Force 
(peak))  Fourads 

.0005 

.26115 

.00)1 

.52632 

.00)15 

.79300 

.00)20 

1.0657 

.003 

1.6354 

.004 

2.2630 

Considering  that  the  magnetic  force  varies  directly 
with  rotor  length  and  eccentricity  from  Figure  (67)  ,  the 
average  eccentricity  from  (a)  to  (b}  is  amaar/2.  The  force 
(F^|  at  one  end  of  tha  rotor  therefore  is* 

rb  *  17.42  (tottotjj)2  a  D  x  x  | 


Satan  GmtaHihs* 
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The  location  of  this  force  is  two  thirds  the  distance 
from  (a)  to  (d)  or  2/3  x  1^/2  »  lm/3. 

The  force  F  is  a  peak  to  peak  value.  Therefore  the 
average  force  =>  1/2  x  F/4  ■  F/8.  The  total  moment  be¬ 
comes 

FI 

Moment  (in. lbs)  =»  F/8  x  2/3  =* 

where  F  »  Force  at  max.  eccentricity  (lbs) 
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